“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1987 


The effect of a turbulent airstream ona 
vertically-launched missile at high angles of attack. 


Roane, Donald P. Jr.; Howard, Richard M. 


http://ndl.handle.net/10945/22249 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
uh D U DLEY research materials and institutional publications created by the NPS community. 
«ili 3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KN OX appointed — and published — scholarly author. 


| I LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








http://www.nps.edu/library 


at lial ct nl ate Lal celle lal li ices 5 sl hol il amar foal pee REPRE ms oo a A a) eel =“ a ys 5 
eH OLAS SAA” SB 1DIE5 LIM Mew od BA AT neath at ere Ce ey ee ee ee ee Se Beta oe iin Led i bn r Qs 
eed eran er itt tei LD) deat eaalehielell ST Pe ne ON ee an Aen aang nai Se aaa fas aba Ya lnateten tea ba eal hy 2 ! be by + yf oe rr Pad te OO a a Sn ttn een tl aaah rE, VON ra RE = SeeeS-ar* 
Fn cnetdadel abel hot a ube ae PO TS eS Cuan UL eee or et (Vn apple ie ore errer rie ee ee eet plas i & Cee Weer ‘ i red H i ol Cv a Pre merexterr Lilt Sarr ect ra eh ene Le ae Po ere ee ery od been en bales 
Se La ee te dl 1 kc Ch Ld cer vena beh deeleeapaline dat ap BP tenh asd Deeeasd alieclalath. ah ele ted pth & Re eee ee Re emeaiied nab ated teent ages or had vs ad ® Ls ee he | rs Ps Cy evs PS ay rere reer re ek ee ee he edna ee ir ea Ba hs oy 
CO Oe nee Te nh cade ete Le te ei oe Pree re ee meet SI rene Crt Lied tAohertc® Saget abn’ eteke% Goes n abaaaitiads Be Ben is ont md ‘ Sd BS ee ey Se a PsP rhea oe earn Ne We Me eee Dk lr eh de ee ee Pr ear telibctontid | oa tiaiitat a: fees 
"] ape eb pare os oar eg oe eps rg. comet WE rubteet bb or PT OT Cte ew eee eee de ei ae Pe Saeed arene ait Sd Sete he fr i or a rs oe em ne vu) mee hs ea caer MTA Vit SD STAT Lo 6 va ee ie er eee tare We PS) oe th ak We he eee hd ae te ae : 
Te ee eee te ee en ee ee ee ed ever ee TB rev ee ey Pry ae | ee ae Pt DE ae tide ae bd = 4 et a Pr Parry 4 en 28 6 ts pre ks, ar es Weer eee mre. Serer Sy ee oa AERIS The ke Cie ele LT BT a ad wevipor cs be bg aoe! 
gentry ae ere rere TTT Pet ee ee eee et Oe et Pn PT et ee aia ee ot Ae Se eer ee aL od La W . 2 ny ri re er ee ec ert re Wee or Per ne ee ae eae ee Sartore err Teh Ter ee eke 2 ea are L dod a asdenpcebe-doe-tonla-tet 
ar et ree eo er eT tt eh eee Oe Ln ayy 0 PY ee eae Pe ed ee Led ra Fad s U a ( 4 . 2 Pa BE Ra ee Rade np ON Ue Met OMAR Tele. He RAN BU Ma astanes ere y tere mee ee ee ee ob pate 
PIC RR Dee tee TIN cote MOTT LO Ser TTY ts Berne Pe tT ee aa aad nee eas Se ele a ee a q : re a a re % « a | a A oe Pe eee ae a Per ee eee ar ro Te One OO ae ee dieteetn Ee Te a aap etna baat 
ee ee te oo ate edad Ree pb site Ra ET Ne eee etek as eae Pa pee et ee eee ee Sn ee ma Ow be . a tte Pees c Weer er rae rer On tare ee) STO Pees Oe Od tert rere ere fee ert he fee i. el le 
Cee ee he le ee ikea © Nt athe singed eo RAMAN UAT -m Sarto IAF et TeV Bea tune! Mat = F a as gra ot Bane e ES A 5 Ye . Ra Se ee ee Perr re kia Peet SSM Pe et ee hd PA Pee Cr CT 00 A ene ence ee de mete eeeee 
ee ee er aT okt ek pe en A el de cd Tal el le ee er T ET Pr he er i ede sa oad shades Lae et of te re " i Pee ihe yo pane ee PL ed es Oe hl Petr re er PT er ene ee en ea Lhe le ae 
ee ee a eee eed a ira al alata Pe ee OS ea td aie en) ee ee al °? s h Pr ee. ee Or Te ee ty er eee Tk td PROT Oy a ie De eh alee Oe ee ee Toa 
er ee ae ad ee oh ele decheAadeas Sib otele a edt hak Poy ree ee SO et SS ok a Py Ce el ee kT er eee ey Bey er Ce or Ce en da) en ee WV ta dk a Vial il Dd) SP ee Cs Pte PT) i bleh de hale eae 
eS Pe ee dE ele lta Prt er eT tere tt ee Ped eee ee eet tee Dy a eed ee ee See oe ee co ad at. rye ahetos Ce. P Prod Pa rT P a fee os ho Sr Cry etme thst ea Bo we bet So a Pe ee a tT ee Le eel) rsegrte Be eth heath eat 
ee ddl a el | tee oad pe eet) Ya Den PP ee ee ee Pees y ) ee re re Le ee, ee lh ale hd bf : . ee ok Ce ey ° , i ia 7 oe aay nena pene Sar eer preemie YT Pt ehh eT PG REL Wi edt deli tell Ce ee Mad den Btn aL ls dali chelate dali caecalieaes 
ani argent ie POE Mente TT ed ae ed an Pearle ata ere ta tee te al ee a aa ae ey Cee Oy ee eet er tees ie | a os . rd OL ed Mr) Pres ere 2. * ler Perr my LY PCO TS TT Ta re nn ta eid tee darth thetati Oe a ee LG Clie tak lied heehee 
Tt ed bl ald ECS ar ep eee ee NS eee fet Ce Co a ee he el Pye Ce eT ee ee ee Per: etre er oe et) oo ey a oa 4 | Ce + p eA a A, ee ere ee wee eM SL Mn Sot wr tn Prva st.) preven ee er ait ere P rer i Oe he kee he ith hs 
a el ed alae sad eked ell ail PY error ee oe nO ee ee ee PO OO ed Ste Dime a or Po eT re = o”~ rT. az ee ret eee te ty eet ee Pree Tee Pe ee Tee TT ht Le dah eal del bled ed Se OO eee thd bid delice adiamdaae 
ata aot am Pe eed Ce Se re = a oe oe , te ® ? das es Pe hee er oe Ct hn hele HOM TAH Se he taMa Tl). Baty WePeeges! 
au rT 


Se Ed dl Oe ee Pe a eT) ee ee BO Or ea ee Pe e 
anaes a BT ae te ee Ba ta a tl Ale OO ee ee Oe ee A Ce ee =f ae a r ee PY he Sr ee he the i te ETS 
Ne lead teh a ee PU er ee re eed ee ee ee ee ee ae a ed Cad ol. el J i % Crt ci RET) parr en te OR ake Re kD 
Te ee beetle a eee head pe ern Peet Otte mr ie ay re ee) ee a a al a Hated t Tea . a Teele tee ie 
rT er ee nT re nen ee ee ee alla Rad dat ee lain tae Raed al eee eel NaN Tee ee ne tcl ' t %eh eres eve ar ne renee WP ee ee Ts alae le 
a ppp nena I Te ele Reh ae haben add oh ok et ce te ai ae al lad acne a a Pty er ers ee ee er ee! Co pi S 9 th Pr rt eee tee eed es ee ene eee eS eee Pere ern ree vem nr etl Po Freres te OO Oe a ae 
ee ee ee Ee Te a a a eal ae eel aa eC TY an On ee A a haiaitdl eee. aes ae D Mtg let! .Ass #- os avr) bd S Ns eee nec na athe ai baie ediral Poe re re veer ete mT Po 
ey el dl al a ahd eet he take Papi peerage NY oe PT aa aa lle tee iseqeng © hateane = oe oat ld a- pdb ot at — oe . Reet ee YY ra t= nb eee oer eerie a Perm term ee vame ere on tt Stet te aL TT ieee aOR A tied 
PT ed eel ee el lee Pek eles ieee le ie ee ee Pian FY ee ee ee ee ee ee Pete ede eee Ge eee ee ee ede tilted ite dll ain a] ro fl 5 ” Sea eer ee aes Ste Om Ace eke Rp Rat Awese Mathys we A OWL HOP ed, pith ite 2 eters el Le 5 4 
Spee eben yr ier ar it See oP OTP E i nee ae el alee Pree Se hee ad ee ee a . Fe Ce -% 4 a ea ease > Pe Wal rT par ner ee pene tar ve mre ea rte LD heeds EPP te ert ener re Pee ee OR CL laa la alate 
SR PT ee dae Ne De eee eet eee Te te id a ee le ea i a aa el pel hen rG able =f om oh aT 0 ee ee og ae A a ew Hl Peres ee Tn Oe de en nec ababend eee hon ee ee a neal be telltale alae tintethet Ol 
CE ne ede ee ny eT Lae pee ee ee ol, ee ee ee ed he Pk So ee er ey e # uO; : . Car) Core ay eo eree eee Oe Ue Sn pr wn ry pares ere Pei Te peperenrty pee me errr es ee ee ebb Pe Se be Ce terete 
a PT Te ee el Mah da ek et heel eat celta ee re eee, oe ee eS ee Le ilies Ce ed a vee» 5 Pe ee eS ee Ce =e o oe Py = Cer rr es Peers eee nr Cree! lnte-ae-evrs ree Peer err ert ek. er tr hb Re Ya ae heat chal eta et ere tl 
a ae Ane a I TIO ed ae a a a a catalase ae alae Od ee ee ee re & Ct tne! co oo : - re ee ere Mere e ta ree Peer ee ee Me nt reer arrmemner nema Store Pre Lett. LL hiannene ard sale are 
a Oe ee ede a ie eo hdd ee ee ee aa ahaha aaa Laced Pool , 5 el hee ¢ J ‘ Fa ty in a Ot Pat a ed ow at. Pon tae eye Pre ar Ten ere rene 7 Oper erent tone atte fr Le aes Ce he et oh tea Re eT eT lied ree 
ee ee ee ee a ad aaah ie ae ee a al el cele erie Ll ater Cd es * . at atm — , f os re Pe Saree ae Cia Se ep re Ge ae EP TR Din dale lane meeps Se ee Ketan letedhs bel adele aS 
apt aieapreeay enat II OO INI TT_ De LADY REED, ET Wat I Me ert te Yd ere OY ee Pan nL st nl aane bean anne aipntind is salah ee Cf 5 Se hata <- en eon Pe ae n Cere re Pry oer Pe Pe Ore ae yO tee Ca ta a del be ee alld Pa PRO Te Te Mea eee Seah ae alia dac iets 
ne ee ee OO ne ee ee Ee ee | eee De Oe i Le te! Pee ere oo Pog eutnat te et ee eee ee Pe, ° sD = . * eo we n we opt eee te ht ir Serie rrr oh ee et ee ee 
Oe ee ae ee ee en ee ee ee ys ee Pe ere) ee al Xe ee ee Lo > rd ? a ed aan! rl i A a heeen Ts eer ee ers a re reeee ere rere YT) Tey ee Ter TTL Ts S Cee re rt hal 
er matt anne Legh Oe ee ee LEIP R I PIO TIE OD Ie Be en NO Oe irre I LYE IE ID I Dee Salad bd eee a jee A ae 4 a ° . ong o Pe Sn ee ce Seehof CP er Oe ale Dandie tha eh dshdalatn titties CT chit vauatedalvniieapussieage tend eateries seta raps teary 
a Te Fe att Ts dad adhiiemteatet ies el db ell eee eee OFF oF ft BOMBED Goh 1 bt atel ote ENS Cory Oe Le Lt ed Prt ete <a a ath FS ed . we se teeeye z 3 re et ere ea ee parr ar eaete are te reennn Oi it Tet tere Per Tay eer le Dp ee Janey CURD AF UD TARY 84g. PE TAoM PE? 
pap a Pe SP ee ed lela tinh rpaceeatiialiad ead we tte = eee pea a soe od ar = ee et ee A Rene ee en PCO een en ee eae ke debaeabeeeter) Pe hii alah daletahad Raila haan eae 
EP eee nen ec a ae en Nata ee a a a OO enn en et ee ae al re 5 Pa Oe a o a ‘i Pee Ve 0 Fe rn Te nett t ddan he eee Yorhlededeebied 
ee Pe en ee ee a Or le ea a ee Or ee ee all nipttaiinnsaatietiende tes | lianas tntaet al ated - - xT an Pee ee Oe ae ee Be ey es pr ee ee Oe le TO a kl eh ee aid Pe eh Dee a Der eta de en hele eit 
Py a ae eee ere eT Py eT ee ee ae 2 Oe oe PP ee a ee ee teed Dae Hi w * “8 ie ru ° ad — ry Pr De) a a RRO. Per ary ort ere Perea, eer rr ee evr Te te a ey napa eae heen = deere eee eee hea eed 8 
Re ee dd ee oe ee ah ee nial ad PO Se ee es ee OT ea] Ne cl ak P olststeaheal ci a) aly eee eT e ww ee pay eet ree ee ears ete TT tea eee Le ee er oy eee ee ee ee eee Pre en hr a Dh hae lanier 
Se Teen tn te a ee re ee Te ee ee ed aa a ee ee eet We eee ee ee eet | des) cd 2 ee = Petre eres ne ornare ie pee eae ery pam rene Ye ee ee OR ce he ea bk henley 
eye er : ee rt err en OR aL lla 








ae! Ce eee ae 8 OD ee ee ee eet hy ek ba 
Be erga rere re wrt mre Srelee Se ry nak Cries ee eC Mae tere iPro Th | 
= et Pre ere er err errr ey Mra te Te oe oe Pr ee oh he dee blaine tae 
See at ee Ee Se Le! Sm LT Le Sn, een ten er ee be tiies 





~ a ee 


bell 4 













PPT Oe ei hh) et th hed) 





eT aT ee eg el a ee ee eee ee ee eT, oe ol eared nae i ee ey PS el a i acd ea acid 4 4 ee fe 5 Pees yar a Serine eer fe nea amr wir PE IMTS Pe et et ee ee ed ad imi 
ee ee ee ee! ee eee ee ee er) Peer ae ere YY CR ee a ee ee pe lat at = q = i re Va oe ere ee Neher On erm err ns ath Ma tT. UL Te Cd | edhe bmeaieh taht hn dnains ee eer ee Ont oe el hice 
Saad eae eee cl eel Leah) Sen oe dd ll Ree ee ee ee ry oe a eed Pe eee ee a Pr ie , el Pe Cha eee ee. bed sy eo a i) en er ee ee. rere rt een re yg we te Py: ee Tre were rol) ane @neyiet ee eS on te Be a Phan dell 
a ee el ee ee ee PO a adel eel alae CE a ed tt ed dad “+ sd ae Sd Fa Pe) vay Seer were cent Sern te ee mee Pere Se a hed a) Piles pa ean ete Ta ah achaislennheatay ete 
Nee ee nen ene ee ae a ae a led ae ae ee ae es ak Nn A lel Ce, el dt Oe Da Ca z dhe oe ee oe ms par re oe a er es Le eee ey eee eaneaee penramnreeey 1Te te rt ete te te Rh bd deta DO Re ea ad en hh Lela aie 
na nian SNe Pee ei OT PY eee Ot le Se Ue nk eal neil wees 2 eee Mere) ewes, Ce eT od = 7 ro Po Ps ” cra rs Tl eee ter ae we emir Pres Ure ee ern en ee ee we ee ian eer EE Tite tn Oh hehe deselect 
a oe el a Pe te ee ee ee Fy ee Oe td Dalal Doh FT lad ' . er ee re Pe ee ee el ee ee) Fed ke Pee he ae A eS ae marta Be re Bs mateo. F=Soee 316 HEFT 
nd na el ek natal et il ont came al ala a nae dint orate? ee ae) Pee Pe Perrys Pet dr reed Ee te Pe Sek Fat? Reg | ee eed ee Cn ’ we Pi AIG c é te S > ~ eh AZASR ee ee oe ee eae ee Os A ks Bt ee ee ae Oe ey Pee Rk ie ae oe eh beak 
ee er en ee ee Ct eee ees || OS i ee Pe Ode ML ee) te a PO ed a De hd A Leet Ce ae arte & &. aS ee ees ee ee ee ee ee ee eer ae ye Eth eh Be aS BFR ee WEY AAA Seets es EDA wes eset 
Eee rian IE te ee ak ee alka eel ea PS ee a a ee ba = ae ld a ah dee rane wf see Fe id © - Bee? ny a. ny Pee an] Pe ee oo ee ET ee PO ee te ed PT eet athe otal beth fen R be tak lathe hele ha aiaeialined Sn le 
- Oe oe ae ee Bel ad Co abe andi hes haha, a a . n er) tee Te oe e a ee Ts a hen pad Pe SD hk eC Te ol Oo eh einen he dada 
ki . ; = : m POT i eng eer tn he endo 
De Le ee Re eh EE Re bh ite tel rt eek dak 


re ee ed ee ee ee ee ee ee eee Le edd St P 
Pe ee oe Pe ee a) a a a ae led 4 ‘rw - e F F ae « - a ‘ad a [en I Pern ee etter Oar ra Ser Ce Vee we aad 
PU er a Ek a aT oe Cry tet ht CL eaten oe Sil hail 
































pe ee ee Nee ed ee ee ee Co ee ee . 
a ale ll teehee ale neath el ete ae Pe ee ee ee ees ee eee er do oe ee ee ee ee (lade LL PT ds ited “ Cn eal " age = oe Oe Ct ae) bh. sam 
eT ed ah a all eae ee Pee oe Pe es ee Le eee Pee ee ee rd ch Se fe.t« i SE ad a ~~) s x . er a nt ra P eee eee re Pen rire Cremer ens ta aes) Ee ee ee a eee aa Se et ik ee ed Se eee blake tah 
ne Na Ee. hel aed ctl al ad thal eth cea lie cele anal OE he el ead Pee re ae oC ee hd AP Ry = ae  ) en, Me Pe , r Re a ot ns - Cee We et et et eee er ee Ce aed Ps On OE eR oe tl hed hal A ee eh eee Lk! Oita 
ee ee eal lL al a alae nl ae ll a ae al On ie ere ne ee be RE ee Ae hh | Ce aon. y CAS AED Oe Pe a i .* * f ee ' Bre? bey Soe ee ne pee heermene an eeer oer ey ee ee ht et ere ee ee aT RO RD hah teed aes el ei 
icp SGD eee ere TT Te ee ee ee ee ah laa ee ee ee Ld PO Gee eee ht ai | Oe ee ee Ce bs re Con ee ao oo Pare me en eGiairereeer wrvamemry saree re to ro aia Pe OR aah Cie eh bed te hte ebb 
eS a een ee ae eee a Da le ee ee a Oe Ne a ol A hl al a ht td . “— ve Yt en Mae dale) ao @ ba we 48 x bs . aw - = pg age A OO nS te ee en eee Tat a dancehall Po ria | Oe EL a ae 
a gy ST Pee a ed ae eS LO aaah aaa Pra Oe | ee Pe or ae oe heed ed eects ad eo nl a ’ ee a) Pre SS iota e = «& OC et et dines al a eeetee-teied TT Toh kh Sa tien Sereatetel dencitinndier tains alee tai sialaiteaaatien 
RT ee ON NN a a Ne el al el hehe head ot ld Ed ad al aed een ee ee ee ee ee CP tL) Oe a) el ye) rtd ona cd Aaa a Some C Tt e a eee 1 ee nie ep ye ee ern ene Pe eee eee ee ld Se ee ee er ee 
ee Te ee a ee ee ee a eee ee ee ee te te) ee Ce PY a) ae i) ” - ne mat a ae wat wa) ee - Pr a pg REE er Pa EP pee tenancies. Seyret pe, ep reener annn meme TT eer E Ye Nae Peer ee 
Od Ld he oe eel PF ee ar ee Oe Oe te ed ee ee ae Me Pe ee le aa Be ed Pe a a ee i bo ow ty ed i Ca te Fk en Pos we Py Fy Pye r ‘as oe ta rq pene “meh Pere Can Note er tte tiie Ca ay Pe eT en a ey ee re er ee denne 
ee ae a ee eal a a A ie cattle dedi L. ° adh baat alia hr ee eid “— Toa Pa Ov ory ran — , A oe er Pe ea ee pea yas eye Te Te ta a UT i hi ine he dha einen a 
r a r D we 2 FY P Cn ee or he Se ea ee ee er MO Se CL aa ee Le Lo 
ee en ea ee Pee hee ee ee a ee ee eh oe Pre ere et ee a ks CY a) eee ee Pe ee | cs ee ae a Foy sd ok J mY ad a me Vee 8H, Tee n = ah om re, PB oe ata eet ka. Pe eee ert eye ae pear ret ee ta ee Te ee ee ee ee 
a pe a re a a et et a a a eT a Ae ee ed Od Sonia * whit ? a oe a etal a ye Td we a eK Po Pe eT et an ar a oe ee Be rae PTS a eo eT le halide ee neteidtt Pe Stel tee) ee he ld hes 
a a ae Pe ee de ee ee Pye Pre See ee a) ee i at oe a nh ee | Fund ¢ 8 enw 4 ee." Fs 1 ae an F, A i a ae ee er A De et <n a emer ertn wer er tir ovine tant) Bene Se ere he bed eh eb Clea! ee tales 
ee Tl alee a ka ciel tl ae ll al laren eta OT ae te ee ea Fede EBs tel 06" a ey et es Cred Pe | rad a Pa wd . ee eo . co ja Oe et er ee ee ee > ee ah Oats “tema ee Brent eee eB lin Nie bate alin lid ee Oe hah tne - Dat ak le ad 
er ene Re a ME a a ak ol aia ie ae oe ee al, eh beh ee Pe ee eee ad Pe re we) oe ee ee as CP ae Py x Oe ree ee ae ee Y) a wa _ - ye spar OR i op eae EE Perreraren pen ree et Paes ee wm, a ate es late te ae 
ee aed sie al atlel athal aheal ed ee ed | ek ee ee oe A ia TS ee al See a ee el eee “ Ce ee we a Ld fo = n ie a ane e P * ort ry - a a. ons ee 7 fe wa ee et ee Ree es PS ek OL Call Be ee te RR ee hee) Pee Rie hie ee Lk ee! 
aaa arn eT re ts ee ea ee et Meee ee Lee Pr ey er eer mrs en ee dal cea) . Ce wig © egies Cd Cad 7 z a " a ny ae oa cies r r . aoe ee es Peer etn ee er Te ee 
pa co e aw ry c °] 4 Ot ae Rhee te en rena ee ee ee eR ae dae be ee it Baten 7 
ee ea LR Ea Oe en OTOL TT er ee fo Tee Nt ot Serie Pe ee TE Oo a a a a A el "+ 0 me v7 *) ~ ee r. Ce ee Pe Regn TS ier erareerry Ora Sr To Nh ete nO eal i dant Lia ahd al lohan te talent 
Saar seria eee wy Pe a ee ee ee | ee De Me ee eM a dd TT i eS Delile iu co eee er ee ee ee el) eT i ch at a re rae aree Cae Venera. a 3 ehrtmate & PReR ee ra tenn rn aCe ak ee eee Se i a sdk oe te Lea eee Peek Th 
el ead a Leal hte dilate tenia ok ed et SR) al a lil lad ae ad i re Se ee en ted. ee Ts ° =t Oy aay a Pee ne ee Ta Tod ta hae os ot PY Me “a a es rere a a a oes A ar ns ‘eo ern Pee ee tn Pee ee ee el Ek Ly PS lia la ee ele lle 
ee ee ee ee ee ee ee ee ee eee ee i td PP a a ee dd Ce al es eo ate ‘ %e Coa es « . a ; ar Pes!» ons ae ay eee = a . mare ae ae pr at bees ar Speier PoP Pe Pan mE La rrr Tt et et eh) ee ee betaine CE te th tl 
Pegg eaten ae yee fener geen te LT I Pee eee reer ye ame eer et Ee 2 ee a le a a adi a Fa ee Cia a or” | Cer eee oes w ae ed 7 “ w « mr a = ee es 7 wag AI Celeb Hse pee Sher beeear een enero Ee a ea ated PP PTS Tena a See Rie 
ea en oe i a ne ee al CO ee dL ala ae ee el a) ce ay oih8 prerin ae ee, ees eee ps ee aay 5 o - ro wr eae ie tar Oy Wi OE a rae ne me ee, Oe, Te Se ih teenie’ eT nt. te eh al tall vr 
Sa decimate ahtilaee tia cated el ah hr cael the a el dD ol ee el aT Nn de a Pert Lad an eee ee er ts ee) Ft nal hl Ot Seer Re ® ay Oy ny ry Paes a PL) a S ar by aera er ary ee) Ryne Sarena yy te er) Ne ee ee ee et Ot tee eh Ne ee LD dale 
dae dee, ol atl aE el cee a ad ele ed do eee ee ee ed ad ee PO Te al et at 8 rH ore re ee eo ded ef ie od te ae o ty tas Ca Pa wor = al a Cee A a | . Pe ey Se aed : pt Care Pere te mete ee ey dn Pe ee Th a De ei ieee lel So behets Se tail +8 
oo Len aiinatomstatintintamped aterm Ce ee ee ee ee ee Oe a oe ae Pe) ee oe ee dl a foe a de lee a Pd _*= yy etd +e 5 t = Pe Be aa fr See eeOETGE ence al srurep she sn \ae re Pheer oe tre Tether rer ty ae pce ee Pre ee a kD a 
a el a ae del ee ee ld ae ee P| Pr ee ee ee ee Care hs oi Fa Py “te er D Oy ' rd aN a oe oe " Ae 7, Oe ees ay ar ve ene 8 é xe Create ee Tre ee oe ee Teer eT ee et rn ee on tt Sen een sre 
ete OP tn dt ee Pe ee Pe en ee Pee ee i ee es ee i ee ee ed = a a) P Ce Fy » - Co en ” ns ry ve Pye ee Ae) fo an ar a we A ae ean aes eas Ws ete Coney et eae ee ee oo ee tie hee oT Ne ce ie ec hee 
a dl et aad Pe eS ee ey ee PP Se ee ee ee Ppa ore PPP ay ere ee et Ore ee me a a | a ae ee ee * rs Fe = - ‘ Pry Pane ~ pow vs ate aa Roe ee eee eh ed oT ae Le EL Ota ho aes cae” Sei teed tal eliate ed 
sarah sie anger lp Tce aeatarnatinn ate teaiattir pine Te atone nee a tinal a itd ae Lo be a? alll ete da tea age he oes a atte chk aces dg iw S j 5 ae ees pepap=anuner ne yeti ers eam renee. ae he Pen ae a aa el a lala dan deal 
wee secmiorel Aan Hemme heated aaah ambukalbriataapemntiarceens pint haltadh wen bitin haGriiiatestndshd Madea ae Retin a eanere = II GJ tn tard eg hi e llth co lla aiach Sad vo LU : oo « Pao Y a 4 Pea eae eres sls nr aclapttee ONE tren Oe Oe MEO ee Od) RR nee fn a in nad oda hicodlirthla endl enna 
ett k diti ie shai el De ell atl se al Pe ete Me te Oe en eee ee dd ed ee oe Oe ee) ee ed oe wt a Oe 7 Pa 7 Py P eee See Poa er par meerrenntan ae OC Ae Tee we ern bn Pa hs Tate tie teeta pnette Care mee ye TH Pee epee tee th 
ee ee ea el ae Md ed ee a Te al a lah ae ee ae ae Sy ee a ee ee a ee Se nd e 1. J A hs we al a . rie . -»y *% q= On eer en Td cae PRP re eee en ee 
Pe te ee Te aa a dN a cd a all ak ad ld eee) Po i bled De a dd ie Ce ek -” Ce oe) i Cy Par , e oe w = ee vs = re eretce “gee Wey re eee) te oe ks ee he te tet Dd ed 
he aE a ele al al leet te at eel Pee ee ee ee Pee Ae ad ' Pd fan wt Pe) e . _—- a ye t% - oP on Pas are elena wes nr Pie es 2 Pee re el Oy ee) fe ee eee ee ded eee ee a 
ee ee ee Pd ee i eed -_ @ CD alee : ¢ a ed , ead = cae > a ro Pe ae etsarte «an yn ey eo re Pr ee ee heh cel oe hemes ee ee ae Dia etal te bell de | 
ee ee ee et ad eed al ool Ee eel la etd PN eld dl i et i a Ma ated 2 ry ‘ Se) A 1 uw. ar ae 5; Poa en cer eee en Yer ree any a Yer oT Lele Oe ee ee ee bila dad tte tal | 
Te einai ah cen a all addicted el adel oil di eT ke oe diel ale ie ici a dl ee dl ee alee od al ee ee Md ras i Py hy =? Cr Pa 00P*le Mme nee age sus or ay my . ne SORES ‘ A & Fs ar eee SS CP ee er sy Sp a ae ee Pees PTT tn oe ee "tee seNene 
Re ee ee ee et ee ed ld ad ed ole EE ie de es fe oe ee a soll PI enme fb Ob a a Ce Sa ~~ al 7 e a . - re, Ay aa ee prepay aseeu- 7 ee ee er ee tee Deo Rel ain del tlh eel | 
Lape See ee en nee een nn een Lapin chbeaeiinatirmiadlntehetiet ii aitatatne I Aelia tat lla ain is oho a ho eel ae is fed Ee (IE ais 8) 8 , ahpint Ls is U x dor ee Oe ee ee ee ie eT er ee en hh ah De led ee eee Sh 
sented batcastlinitineadll nathaghadeenradh etntindanntiadpodrghedinats Weaphatiale aia gha ea fisiid, satiate Staats a ee PM DN i eS aol ache cal el chad hte bah iia! otal nail Loaded Fd 4s ra er 2sh on” er aT ae ~ ee ere en? Sey IT el nlite nae aie a at een lan be a 
sehesiedoerall ped naptaacetiadedl ceded paiegli ah od ebeenteh Aten aasetteh Sake eee ee se ab eee lana lll ivadiiand al i = Mle a Oe ay SPL ae ad a * Ps = e ey Or Oe a) Wer we SC Oy WO ddeedlaaetiee aS En aie 
eects. shat tapn retaelaln eatin ie meneae, 2h hentia) tog $field ail ica dd HO >a beep idegt i ai =a S en Nees a! An bade ee] a a a Pare nea rent eptar aera et Py RR OB Perey ie Tae TO ete alah ete she ie UB er etereen sr eememsaewvatrnsy 
_Aeiatall eed tan vn tial do me gg PORE Pre gener ery ft ed , ae ee ee ed Pe oe rd re a aa ee ee Te te ' a Sd i ae) ° a 5 -, eer A rere ae as = ican FORGES a ery ay ene en ey ee ere ent oot] whe neh Bg | Te en bat MR e eeeen tens 
Ce ee el ee ed ee a oe ee i lee ar et a Pe ee ee ee ee ee Ta cr ie ae es nd oe a - ry P , . abe fd = : oi oe aaat “ores ma en et eon 
- be A Pn 5 A * it ae eal be ae ah al zem DO ad Dey Cal Oa a th lot athe tet le Dee athhin deiateg! sh F 
spbmanntretns! cen giecies ee I a alll tee ee el tl ed! ell ed il ee ol ied a ed Can! ant “et es oe e Coy eg o , - . A he ees eo Me er es _ = pore err Pe ea Tre) ed eT ere he ee el ok at 
Sie det ee a ee ee =o a es eee eee ee od Pra Pl ae * od a LIL > wet a Deed . ay ‘ % fi wi 7 Sok a Te ee ee ae ar reagent tne ee Se a Fd an hie el alee ati ends 
=a gots beeiales Mill ell acento hth atielleed ad ctaha ed ia el Bnd ae ee al - ol etd rea ass . ae ee tt “ 4 ” Cy Po « ae Pr ‘~ ax Seon er wan ennae Pee yar A er ee ee ee ee Ts eh ain tLe 
Lipattiipel hae Miah rita aes ty thenficlheehiadind pO hia ei po ee ae ee Pe ee nL Ne ol oe, it ae centile ted sad sh c.e a “ , t- F oe Parr. ba % — 7 ee eS rt ee Oe My Be SE A tedetadambechas hears tlhe 
Dl nd ee ee ee oe ee Dl oo le ee ed a ed ee me To Pe ee 2 “ww v , oo w- ny a a hee Se ee ~~ i Se | a ieee er re er | ee pe ee OE a Lhe aa ad aed 
any an a Pao Oi ee eens nan Pe ee ee ee al es ll eal oe a ad ed we were) ’ Sa a C mi | 7 > he = A a Ba pepe ot DT at en Dede ated Tes Telia eh ala LeD ona eee slieahia Ca = pe ald 
ee a i ed de td ee et Pe las Le el C2 rds ys Ca | ts Th * Ped cere 5 caer ea ei ana er Deter Sor rere Tee EE LO alee Sane er Ee ne ett ttiaibabeeedts 
a kT al tlie teeta el ln ed ee ed ee ee er ee ee “Bae sollte olla oh at b nti St rot a 5 pw a a oo 5 Pi oe. ae OS Cn ea Pee ye are i er Cer ie ve ee ee inte eee ee et TT ele dedi abi 
Ch ele alien ee Dd ee Od dd ee a eT od Pe -—- « = wero in. CY led ad eh alle od fe Woe aun pe See ea sau aga a ae Aa a ae a ean its a i aha fa pr Ree er RS 
ed ee ee ee a ee ed -@ ee ee ey 2 rn les 5 SE ie Mean el Cio.) Ee ‘ ies ite why PS rd Cd a8 8 ry “gt o- > eh “ees eae ape oT es Pir) ee a ee ar ee att i al OT ed ell Seine Meal catatad 
tattered oe hd ee eee Pe ee ee ee ee ye Pree yk ety ed a ed ay ar Pre) - Pe oo rs ry - r id — ee ro et ae xo. a ore ace were ter ee et ney ee ed sa." ee ee ee 
A al Tt ee eo eel ade el ol oak deel ol a ew Oe de eee ee lad Stall . or one . . a en = a ed ] -¥ Par — x a ; e 5 Baalgers mii ry oe om ; wotn te ee = = peels a a Aap aa RE EER eee nae ee agree 
dad oan te ot ee re re a ee ee ee a de) ee ee i ee er deal . Pa = ny od : ae 4 - a ae Pe + balers eee ie ace feet AE pes er a een oe ee ee 
: Saat Shae ne = Pe ee ee Te er eS Pe ee el Oe) dk ok eee ee ee ee Peer mete Ae ac ary ee ny re Ps e fi ps -a ean a arg ee eee ee eee oe a ee ee Om ee eC Le heeds 
[dilate eas Oe ed 68 Plo - et eal ee. ae oe) Py eer oe Pe =. « rary - eo! « oo eee = ee Fs Ps oe < FS a ban ince came Se aR rm my ert Pers rr ene er wernt tad ad 
Geant titties a at Tk dello Di ott tn el el iia NT et eel tal oe dt ee # Pe ore ili 4 ro od > of CO e Pa 1 e i me 86 SP othe “ ood ts ae ee Pd ey Sd ~ te ed endhema tenn 6 t itadiadee adil 
soit -salh ee a ol Ce eld al lle ola atta ad we oY ny a Pd ad ee Oe Ad ed + sa . ny = OI um 5 aN 7 ages ont oe nt = ee 7 ah re, ae ~~ ee ee a) a Re rk et) 
Con titel tanta tare betel td ol oe 2 i red Se a ee ee od do ” Or ad Cth od ee ad ” Py econ 7 eine eee Sha open tera Par PP aaa er OMT 4 aoe Ce aad 
hesitant Se ial a ee ee ee) ee ee re ot et fn wee -_~ = = a ae os aed Ca ~ o™“ ar 4 = A eens pentane te See ee Ye a Oo taal re) Se 
i en a cel ee ee ee ee ee ee) pane £ Per ee | Pe ee e - Se ee a . ad a] rr) @ rar - Ps e C baal * = ° Ce Cot ‘eo babel ~~ = 5 
a s a e - F ig cs chic D ON Sed oe ie ts « “fe ne TT” bm at ne dol ei Se a on rT a (a Sead tn fin_ ind haa 
i oie ass Oe ahaa Rate -_ = ett) Ct aed od i! -- ad aad stil ak aasligs ‘ Cd ° cg C _ . h Py PSS ao ae reer = om a ease ee | ol ee = ST ne hela 
en ee eT ee A ee ee ae ee mao a een Pe ee a on ary Pa A 8 ms ow a Pe Ss Po a ah pe ee ee a SBOE aA Re Been tate a Bee 
et Tonal Pee ee ee eo 2 Pe ed oa --« Oe Pe es ee et a ot -* a Td Fd a ld a Py eer) eS a? 9 ° a Sas ae sie F a Sete a “ards sll gall Seo Pe ee apr ae rea nT hi tad. 
ee ed ee al -PF 1 Ps > ~ we or] ry oF bd at or re Po - * Ps ie . 5 ‘e P adie St a ats Co ks tial ila hal phe we one hy ils Paap eee ee oe 
ee Pet me FO em ene ete wr oe a aed * ee ae a ar roe. nS ow ys - = a ees 5 ae be eda = e jl $ roc ibd beara basdlty ada earGe ar wre en OMT Par nr 
Fame ater tte we oe ee ee ee a Pe Oe PO a Pee Py) er) rd Pare ee - «= - ry ny ’ te e a a a oil be Mi edi dhs 79 aes Seba tpl iret eS Saar he he eee teen 
° ha cd a - . Bees = a . é “ers art -- 
Steel atte tt aaah diel ae) a « Fathead Ce Lee ted ee eel - a ewe A ato HF ed ed Lo eho tae = - ¢ we Ca ~ em aes : 9 a - eee a a. a a “f ages. Peete me ory raat OD be nathan Reta Ld ae Sod 
slatted al “oe Cal ea ol dT al er) n ed tl Ce a ee ed ote = LE) ra co] 4: r - a a ae i= 5 ee Pea ass . en “pe Sala (Ss ds bed Some oy eae Oe ee 
- o ~ = = me a Syn ‘ a 2 - oe Es 5 = Se . ° - 
mien 2 a — ca om rere = Ot ld - od * Lot Ces - - CY ° 4 - ~ of C7 ad ere * = bd pam a = = eye pe Se are ee ee er = Ce OE ek Leal lied See tiieted le 
-_— ec ows mo i i ee) Pe - we a ar a - ee el a ed On - oa “ - a rs _ x x a = ie a ea a oa ae aia eo 
we -_ Oe ad = oy -— & "2 Dd ed a - ed wae - ce - od ee ned ra eT , = e - v ed - Lt 4 . Ot w as ry = Por e ay Se lens Ald a lt Ee at ee tad De I eh Pe no ln tale tea ead = ane 
= a td Olid Ce ae) a 22 | ee use wo ~~ ye © a - fa a a oa a a - ry ~ cS =e « x re apace ie ma Fos SE el = = Ce ee e oe Sh 
CE died a) a) a nd Po = 7 o- 2 “Ch oe ae on a = ete = F e a - PY w - 2 a - & » a ‘ ra is CL Ls Peas eacleer me et ack Pe 
Lo eal oa - - = - = = eye) Py - = aon Py ” e re oy - Fa ae ad ad 4 me 2 6 Sh pee | ry a be he Dl eed 
lat = « ad Pd eed Fa ad ° oa - CY aa) Ped ry ny oo « vr) ‘ r - Bs Caneel os a ae ars es = Pah Ges AE ates Py Ae 
ned ai Selig e Ps i bt — ™ Stara = = . bes ig tz aes ha Li pl P ae ° ra een =<. . i a | = a ee ~~ foe te ee | te - ed ee es ee ae ts 
olla -_ -=- - “ 2 owe o - ary er ar wee a ee) = an) 4 7 Pee a z Pa BY ere “ a _ ee ee oe TT a el le 
eee 2. - — am coe ed ad a PI ee =" a «i ed - -? = rt; a an) *%~ ’ y o = es a od e a Pr) yo ye ee be ie) he tel tl oe ee | inne taht » ve 
a Sener] e Pa bee tae aa oa Poe a aa - a Cece 4 e a i ih : lia 7 . : Pa - = Peg tis Pn ee ns SE ae ea aa = © 
= - S a Ly - oT - oa - ms 
en ene = ee ~ a aa! af a is bs an ae sf ~ B Py Pod os cot ry FP er ry eT) ee ee ee ee ~ 
5 eee - . - - ae =. = = S - . . = : se - x Sey ee itg ee ee - = apres 3 ae il eiaeaer poe Peete eT a an eae ge 
- Se Le os er r oo Ca ae - - 7 - = re B, ‘ es pe teas EEE aa a pes ee a ae 
- ° ee dS -s -« Po a 4 a a] . reer ae. Ort) ua - “ oo - - a = Ce a < ee eas Pe Serio oe - ee 
Le eet aun a a. S 2 * Sf - Sy Co - *£ at  - Ca eee == =e er * a ea 
7 oe = oe acl a a Ast = api a Me a =e ame Sf e a 3?) a " o Py aoe . Pd - bd ee eee ee en a a ee re a - “2. 
ee ee - Pere - he . a P tL od = 7 = B = rs rent Pe a mA Ca ps rs es <= ar Fie Te 
- - - - = « - - e i 5 Ps acum 
ie = eg adil = dL Pact ead 2 ny ma pak -” ae p 9 rs rs a S = s = oo s a A ct —e oe ear a teh Lot ate Me Dalal . na Sane 
Fe E is ze h.) pied = 5 a o= S P a Pa Sere Te od ot 
~ en eae w - a ° Yh a ead Cel aS . el aa a a! = rw ad a id 1@ > . - pr es 2 ia) iat - ° = ee = - paid = * = Pye ages ap pe ee ee 
- od ay * "] oar Ca - a ny Pa - a ¢ e- « s bs - - _ a 4 
ee PR ea Pee 3 zi pa pe ee api Sas ae a i = Cd bat ca Co - ie = cd Cy - Pen oe e Re es >) aol _— ort a ee ee i 
- S oy ad - wen we -- = . FS io - s s 5 * ~~ = = Cad 
e ae a ee Pag AF ae Ps a : mee poe ys : ang p = F Fe e s 4 a F os = 4 a a a a == — a wen Sn! > r a Se hey gleams ae a pare Ee _ eae 
a ee =, Pee Pe ee 2 o Pea ae -y - Pr ry a . ny “ cs na bes Le * = ee ee ae 2 ‘ane pre Tees eat ee ie 
e - ayer See ar an - - ae ry Pa baled ri ny . ae - e os . aT Cr heat a ks al C 
at bs = 7 cd e- o ~ - ~ va a ey | a me 8 ae Cia le i tians a Sete) a 
* - - ° oe ees o = a Oy a fs “oe A a aor oe a ey 
Poet s “we = ° “ x a: re ra as - . ~ ° ee oa Ss a yi G ” = aint,’ haa py od - oi we ae z Ars - ae 
oa —- a Se en ae ~ s E i i Sars en a a Pa Ps « ” a) +: ee di i a be G s a = ee ee s,s +. ae 
= i is Capt - =# i = bad « = aad iad a = = e M MS a <i cad ¥ ‘isitieg Z Me nee be si J Oe ee ee ee ee ao 7 
ltd - eld =e orem oe = - ~ =e = a ars - cars ny ry = r = ee a rd o Pas ~ a a a i 2 - see cS a a Learns 
a Ces Cg ord Pty - - Po oa a i = ra ; - e wo? . - ta - -¥ and ¥ de ey - x A as * =e 2 ae 
Ce OE ta ine ac) a. = ed dad = o - . a - od ° bet a = cI © w* Prd * = Cis - a 
@ ~! Sy See Pe eS oe = Gis eter “ae Pre Pe f ” ~ e Ee ° - ee er =e ee es - + ae - - oe Pe at tds SoC a, 
w--e on a ars Pe ae Ps oo - - er pera is A . 2G a) ee ee n 
Oe rd a = 5 - me ioe ie S . ry 5 - Ld ey = se ci p= ate m Pe a rn ee Co pone ees = a ee ie we + eo 
me Gas » a bo ~ a oe ee eee ~onn ot cas 
Pe ee eS oe / - é % a a. « - 
adel rs ees ae = bd » = ° -~_"= oe oP ee Fa BP CS - ate et 8 eet o ee) Ce he ee ed - Ca 
z cal ns Lat al i! pa dik © ben Lo - ad 7 - , z F. is = eed sf CJ e ee ed Se ale 
ca - td a en eet 2 Sec wwe - Pye a) - e- . * - = a = = eee la eG! Pi Ala Midis M DL Ree ere eCy ara A oer ees - -= = 
Ce Na id er ee ed w= = a] a ae Se ee I ba as A ry ~ er - + Nee eee or) rr er asda a ae cir aye eas o id ay aan 
a dae a a - od - Poe es ‘a 2 = Py Spear Fs . - = P P -~ a Py A aor? . Cho tied — ae eas ri 
= ae sh coke bad Cad cd =: 2 or . = Ss - 7 ae a a . . oes FA ba Ci - fuer ry a - a oe ns Ces é - Pa ate ae me si 
-— « ad a -« a 2 Pa a ed ra a “ 7 ce se = a - LJ os - a = phos pate a] eS ed - “~ 
a - -w -- a) Pa - Fs = as - a) ry i « a = a eae * eo bs he plod oD st a8 _ GA Fie Pr 
ro oe = 4 es =. eee = = « e ar - oe Fy = o - Lo] ry eae are oo - 3 
co Ps - a - = oe 4 = a - - FP = - - a -~ at = ean i ee ated be a 
- ~ eS os - i Py = 2 Ca’ - - on a 
. a oa ee ~ - Ps = * i ey oa - rr . ~ 
Ca Sd - - = - ry . z oar - Par 
aad neo od a . - . - a a ° ~ a aan » he 2 Gat oi 
- oF Sn - aa e es Py aS er ' zs : F a < 5 - - é, =~ ry eI ray fs tae <= reer aed rh a - ox. Co 
- cal neath al - 
- ogo) - -— «- o- a = - ry Par) - By - e - ° a acre 
oe ro o- e = = * ry an’ i hs bat Sail pa = ww a a 
rs = z es es é . -_- = - a 
~~ - ad -— . ad oy - - e ory a7 T A a —s -_ 
al - oT . a rs = me ° = e - 5 Sa 
dl - Py -“ . Ps er A - = FS 7 E oe a a ay 
od rs 5 = e - Ps - - - - as 
- = a a a oa ded - - ey ca < ny Ca a a D 
« - . ny — = a es = ~~ 
a = ra ow - - -_-. e - co a Es 
LJ - - Ld - 7 bd 
« - 
x is 7 s eg ae ° es " 7 rs ae 
Cel Cl * - Ld = cf Ca sl bn rd = = is 
- = - - - = ~ ” on) - A ro CT - - =) i= 
r S aa. = a aie < 2 = = A = = =<. 
co - a a e - - - - Ps - - bed 
na - ra = - 
- - is = id e aS we == 
ke ° - . - « ba - Oa o b 
— - - ‘re . Pry ~ v. ary 
- - - . a A 5 - A . - = 
a Py ‘ Pr e - es e 
. ° - . - =e - . a a oa 
rs - = . 
- es ry . co 
~ 
Ud 
a 
Py - = 
° at pe - 
; = = Te Lenebmnaiin die Atenateaneell bl add 
altel cs 2 Si ad 5 eee ee ied - - 


= Pepe met - ie mn 











NAVAL POSTGRADUATE SCHOOL 
Monterey , California 








THESIS 


THE EFFECT OF A TURBULENT AIRSTREAM 
ON A VERTICALLY-LAUNCHED MISSILE 
AT HIGH ANGLES OF ATTACK 
by 


Donald P. Roane, Jr. 


December 1987 


Thesis Advisor Richard M. Howard 





Approved for public release; distribution is unlimited. 





ONELASSIP IED 
SECURITY CLASSIFICATION OF THIS PAGE 













REPORT DOCUMENTATION PAGE 


1p RESTRICTIVE MARKINGS 
N/A 


3 DISTRIBUTION /AVAILABILITY OF REPORT 


1a. REPORT SECURITY CLASSIFICATION 
BINEGEASSIFIED 


2a. SECURITY CLASSIFICATION AUTHORITY 
N/A 
2b DECLASSIFICATION / OOWNGRADING SCHEOULE 
| N/A 


4 PERFORMING ORGANIZATION REPORT NUMBER(S) 5 MONITORING ORGANIZATION REPORT NUMBER(S) 








Approved for public release; 
Da St mroueiton unlimited . 


'6a. NAME OF PERFORMING ORGANIZATION 6b OFFICE SYMBOL 7a NAME OF MONITORING ORGANIZATION 
(If applicable) 
Naval Postgraduate School | Code 33 Va Oorornaditiatemsc loo. 


6c. ADORESS (City, State, and ZIP Code} 7b. ADDRESS (City, State, and ZIP Code) 
Menterey, California 93943-5100 Weionterey, California 93943-5100 


| Ba. NAME OF FUNDING / SPONSORING 8b OFFICE SYMBOL 9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (if applicable) 





10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT TASK WORK UNIT 
ELEMENT NO. NO NO ACCESSION NO 


Jee oO Orn Leen rT ARS hREAM “ONMA 
VERTLGALLYeLAUNGCHEDsMISSLLE AT HIGH-.ANGLES. OF ATTACK 








8c. ADORESS (City, State, and ZIP Code) 





12. PERSONAL AUTHOR(S 
See e Roane, Donald P. (dr. 


faa TYPE OF REPORT 13b TIME Co 14. DATE OF REPORT (Year, Month, Day) [15 PAGE COUNT 
Master's Thesis 0 —_———— 1987 December v2 


16. SUPPLEMENTARY NOTATION 


(11. TITLE (include Security Classification) 


COSATI 18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
Vertical Launch, Missile, Turbulence, 


HigneAngle of Attack, Vortices 


19. ABSTRACT (Continue on reverse if necessary and identify by block number) 


The purpose of this thesis is to initiate an investigation of the 
response of a shipboard surface-to-air missile launched vertically into 
a region of turbulence at high angles of attack. A review was conducted 
on the effects of asymmetric vortices, turbulent flowfields, and the 
marine atmospheric environment on a_eslender body of revolution. 
Turbulence mapping of a wind tunnel with installed turbulence-generating 
grids was conducted using hot-wire anemometry. The resultant turbulence 
intensities and length scales were analyzed as a function of the 
downstream distance and the grid mesh-width/bar-diameter ratio. 
Turbulence intensity was found to decrease, while the length scale 
miemeased, with increasing distance from the generating grid. Both the 
turbulence intersity and length scale increased with an increase in the 
grid bar diameter. A generic, vertically-launched surface-to-air 


20. DISTRIBUTION / AVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION 
KXUNCLASSIFIEO/UNLIMITEO [J SAME AS RPT = Clotic users } UNCLASSIFIED 

22a. NAME OF RESPONSIBLE INDIVIDUAL ~ 1226 TELEPHONE (Include Area Code) | 22c OFFICE SYMBOL 

| mrenara M. Howard (408) -646-2870 Code 67Ho 


_ DD FORM 1473, 84 mar 83 APR edition may be used until exhausted SECURITY CLASSIFICATION OF THIS PAGE 
All other editions are obsolete 








2 U.S. Government Printing Office: 1986—606-243 


1 UNCLASSIFIED 


UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE 





teem. Lo. CoOnenilce 


(VLSAM) model was developed and constructed for further 


missile 
moment, pressure and flow 


wind tunnel testing, to include force, 
visualization studies. 


UNCLASSIFIES 


2 
SECURITY CLASSIFICATION OF THIS PAGE 





Approved for public release; distribution 1s unlinuted. 


The Effect of a Turbulent Airstream 
ona Vertically-Launched Missile 
at High Angles of Attack 


by 


Donald P. Roane, Jr. 
Lieutenant. mited States Navy 
B.S., United States Naval Academy, 1981 


Submutted in partial fulfillment of the 
requirements for the degree of 


Vio be ne OlmroCle NCE lyre NGUNELRING SCIENCE 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1987 


A) 


ABS T Reveal 


The purpose of this thesis is to initiate an investigation of the response of a 
shipboard surface-to-air misstie launched vertically into a region of turbulence at high 
angles of attack. A review was conducted on the effects of asymmetric vortices, 
turbulent flowfields. and the marine atmospheric environment on a slender body of 
revolution. Turbulence mapping of a wind tunnel with installed turbulence-generating 
grids was conducted using hot-wire anemometry. The resultant turbulence intensities 
and length scales were analyzed as a function of the downstream distance and the grid 
mesh-width oar-diameter ratio. Turbulence intensity was found to decrease, while the 
length scale increased, with increasing distance from the generating grid. Both the 
turbulence intensity and length scale increased with an increase in the grid bar 
diameter. A generic, vertically-launched surface-to-air missile (VLSAM) model was 
developed and constructed for further wind tunnel testing, to include force, moment, 


pressure and flow visualization studies. 
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l. INTRODUCTION 


A. BACKGROUND 

The recent introduction of a vertical launcher capability for ships carrying 
surface-to-air missiles represents a major advancement in weapon system reliablility 
and flexibility. Conventional trainable launcher systems are extremely vulnerable to 
mechanical failure: in such instances, a fully loaded missile magazine is rendered 
useless. However, each vertically-launched missile is fired from its own storage 
container, cffectively increasing the redundancy of the launcher systems. Additionally, 
a vertical launcher svstem has unimpeded access to any of the weapons in the 
magazine; conventional launchers require cycling of the weapons until the desired 
round is in position for launch. 

The development of these new launcher systems, however, also creates problems 
not previously encountered with conventionally trained launchers. Launching a missile 
vertically into the open ocean environment exposes the missile, while still at low launch 
velocity, to potentionally significant crosswinds. The result is a missile flving at 
relatively low velocities at an effective high angle of attack. Under such conditions. 
largely unpredictable side forces can be generated, which pose a serious threat to flight 
stability. 

Additionally, the airstream close to the ocean surface and over the deck of the 
launching platform may be relatively turbulent. In determining the effect of turbulence 
on boundary layer development, both the turbulence intensity and turbulence length 
scale must be considered. In the marine atmospheric boundary laver, the turbulence 
intensity may be relatively high; values as high as 12 to 17 percent are possible. 
Additionally, large scale turbulence (of a scale comparable to the dimension of the 
affected body) is present in marine environment. Large-scale turbulence has been 
shown to be a primary cause for the rolling, pitching and yawing motion of bodies 
subjected to such an environment. On the other hand, the degree to which high 
intensity small-scale turbulence is present in the atmosphere, and the consequent effect 
it has on missile boundary laver development, is largely unknown. [Ref. 1: pp. 3, 34] 

The characteristics of out-of-plane forces and moments on unyawed slender 


bodies of revolution at high angles of attack, caused largely by asymmetric vortex 


shedding, have been investigated for over three decades (see Reference 2). Much of the 
research to date has attempted to model or predict the flow about such bodies and to 
examine the effect of design changes on the observed flow. Experimentally, several 
methods have been employed to analvze both the cause and effect of these asynimetric 
vortices. These methods include direct measurement of forces and moments, pressure 
measurements, and flow visualization. Force moment measurements permit an analysis 
of the magnitude of the disturbances, while pressure distributions and flow 
visualization enable the cbserver to study the location and development of the 
Be munetric VOrtices. 

The goal of this thesis is to research and experimentally study the aerodynamic 
characteristics of a generic, verticallv-launched surface-to-air missile (VLSAM) in a 
turbulent airstream at high angles of attack. Although extensive research and 
experimental work has been conducted by the author in the development of this thesis, 
it is not a conclusive study of the problem outlined above. Further experimental work 
will be conducted at the Naval Postgraduate School in a continuing effort to 


understand the effects of turbulence on slender bodies at high angles of attack. 


1. Asymmetric Vortex Theory 

The flow about a pitched, unyawed blunt-nosed slender cylindrical body can 
be characterized in part as being similar to the flow behind a two-dimensional cylinder 
at an angle of attack to the flowfield. As shown in Figure 1.1, a separation “bubble” 
Seists at the nose of the body. Further down the length of the cylinder. two 
symimetrically disposed vortices begin to form on the lee side. These vortices are fed by 
the vortex sheets of boundary laver fluid which have separated from the body. Further 
along the length of the body these vortices alternate in their separation from the 
surface. The result is a side force on the body, relatively small in magnitude, which 
appears as a consequence of the steady asymmetric vortices. The description above 
assumed a blunt-nosed cylinder. In this instance, a nose-induced flow separation will 
occur which prevents the formation of vortices in the area inimediatelv aft of the nose. 
feel. 3: pp. 751-752] 

With the addition of a slender ogive nose, the blunt-nose induced separation 
will not occur. Consequently, large asymmetric vortices now appear at the nose of the 
body. The aft portion of the cylinder will still produce a side force as described above, 
but the magnitude of that produced by the nose is a much greater proportion of the 
overall side force. The net effect is that of a slender body with vortex formation now 


occurring along the entire length of the body, as shown in Figure 1.2. [Ref. 4: p. 99} 
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Figure 1.{ Vortex Flow About an Lnvawed Blunt-Nosed @yimieem 


The point at “hich the transition from symmetric to asvmmetriowyoncs 
shedding begins is dependent on severai factors, notably the body length. the nose come 
apex angle, the angle of attack. the flow Mach number, and the Reynolds numipes 
(The effects of the Mach and the Revnolds numberewall be covered in the foley 
section.) The transition can actually be divided into four separate regions, as shown in 
Figure 1.35. Each of these regions can be described by the angie of attack Game 
characterized by the type of flow present (see Reference 5 pp. 22-23). Regime Iles 
@ < 5°) represents a classical potential flowfield with no flow separation. Regime II 
(5° < a < 20°) shows flow separation in the form of two symmetric vortices Grtiemee 
side of the body: no side forces or vawing moments are induced. In Regime II (20° < 
a < 60°), the vortices are now shed asymmetrically. This effect produces botmere 
forces and vawing moments. Regime IV (60° < a@ < 90°) resembles wake-like flow in 
the region behind a two-dirnensional cylinder. (The angles given are approximate and 


do not necessarily represent the actual points of transition.) 
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Figure 1.2 Vortex Flow About an Unvawed Slender-Nose Cvlinder. 


On a slender cylinder with a relatively fine tapered nose in a subsonic, low 
turbulence (T,, less than one percent) flow, asymmetric vortex shedding usualiv appears 
aoan angie Oy when the angie of attack @ is approximately twice the nose semi-vertex 
angle On As the angle of attack ts increased above the angle Q.. the region of 
symmetrical wake decreases and the asvmmetric flow pattern predominates. The result 
of this increasing tendency toward asymmetric flow is the appearance of out-of-plane 
forces and moments. Steady or unsteady asvmmetric vortex shedding, or “flow 
switching,» can occur as the angle of attack increases. Steady asymmetric vortex 
shedding produces a side force Y which increases to some maximum with increasing 
angie of attack. The magnitude of the side force then decreases as the angle of attack 
is increased. As the angle of attack is increased still further, steady vortex shedding 
gives way to high frequency unsteady shedding, and the side forces disappear. 
Beet. 7: pp. 1-5] 

Detailed descriptions of specific experimental results were intentionally 
neglected in the above descriptions. The behavior of the asymmetric vortices is well 


documented for a vast number of models and shapes, yet their cause is still not 
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Figure 1.3 Vortex Shedding on a Body at High Angles of Attack 
(sce Referemec a: 


positively established. Most theories suggest a relationship between the flow 
separation along the body and an inviscid interaction of the vortices {Ref 8) poms 
A prominent theorv, put forth by Keener and Chapman, is that minute variations or 
disturbances at the nose apex lead to the formation of asymmetric vortices [{Ref. 7: pp. 
2-3}. 

Several important generalizations can be made regarding the effect of 
configurational changes on experimental results. The addition of a long cylindrical 
afterbodv, as well as decreasing the nose apex angle, will result in a decrease in the 
asymmetric vortex onset angle 0, [Ref. 7: p. 5]. The placement of long strakes along 
the missile body and nose was found to reduce the forces and moments generated by 
asymmetric vortices [Ref. 6: pp. 5-7]. The use of small nose bluntness, nose booms and 


boundary layer trips has also proven effective in reducing asymmetric vortex-induced 
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loads [Ref. 9: p. 552]. Other methods are available for reducing the formation of 
asymmetric vortices. Those listed above, however, are relatively simple to implement, 
both in design and experimental work. While these methods have generally resulted in 
a reduction of vortex generated side forces, they do not necessarily work in 
combination, nor do they alwavs give consistent results, even under simular test 
conditions. Additionally, the effect of the changes may adversely affect some other 


areas of concern, such as lift, drag, or stability. 


B. TURBULENCE 

The term turbulence is best characterized bv the presence of random, short- 
duration, small-scale variations in a flowfield with a given mean velocity. In 
calculating the effect of turbulence on a body in the flowfield, a comparison must be 
made between the scale of the body and that of the turbulence. In addition, 


consideration must be given to the energy contained in the turbulent flowfield. 


1. Turbulence Intensity and Length Scale 
The turbulence intensity, T,,, is defined as the ratio of the root-mean-square 
(rms) value of the fluctuating streamwise velocity component to the mean velocity 


component ina flowfield. Mathematically presented as Equation 1.1, 


fea ul (eqn 1.1) 
the turbulence intensity gives a measure of the relative magnitude of the velocity 
fluctuations. The turbulence intensity also yields an indication of the kinetic energy of 
PvemuureUlence, Via the term’ u in Equation |.1. Therefore, high turbulence intensity 
correlates to high kinetic energy in the velocity fluctuations, i.e. more “turbulent” flow. 

The length scale of the turbulence is a measure of the dimension of the 
velocity fluctuation. The effect of the length scale on a body is a function of the 
relative dimensions of the turbulent eddies and the body itself. The length scale, 
despite the inference of a purely spatial measure, also represents the temporal character 
of the turbulent fluctuations. An increase in the spatial length of the turbulence 
corresponds to an increase in the time the body is exposed to the fluctuation. 

The dimension of the length scale relative to a body can be regarded in three 


road Categories; larger than the body, equal to the body, and smaller than the body. 


Pi 


For turbulent eddies of a scale much larger than the body, the effect is similar to that 
experienced in the non-turbulent flowfield, where the deviations in the speed and 
direction would be of long duration. Thus, the motion of the body would be affected 
in the same manner as one in which the flowfield had changed for a long period of 
time. [In contrast, the unwanted rolling, pitching, and vawing motion of the body is 
caused primarily by turbulent fluctuations of a scale comparable to the dimensions of 
that bodv [Ref. I: p. 3]. Both of these instances are sinular in that the turbulence 
length scales impact the motion of the whole body. 

Length scales of a dimension much smaller than the body, however, have a 
distinctly different effect. These length scales have a magnitude comparable to the 
thickness of the boundary laver on the surface of the body. Therefore, the eifecumer 
small-scale turbulence 1s on the flow over the surface of the body. The development of 
a boundary layer and flow separation over the body can thus be greatly affected bv 
such relatively small-scale turbulence. 

The length scale of a given velocity fluctuation. subject to interactionmaied 
other fluctuations in the turbulent flowfield, is not constant. As described by 
Bradshaw (see Reference 10), the effect of vortex stretching tends to decreasemaie 
length scale of turbulent eddies. This “cascade” effect is caused by a strain in one 
direction affecting the orthogonal components due to the conservation of angular 
momentum. For example, consider a vortex rotating in an x-v plane. Strain applied in 
the z direction causes an increase in the velocity and a decrease in the length scates of 
both the x and y components. The increase in the x and y velocity components then 
similarly affects the velocity and length scales of the y-z and x-z components, 
respectively. The “cascade” continues to decrease the scale of the eddies) Wntiiaae 
energy contained therein is finally dissipated due to viscosity. [Ref. 10: pp. 14-15] 

The cascade effect therefore produces, from a single source, eddies of various 
length scales. At any point in time, turbulence may be characterized bv a 
representative, or “average”, length scale. For turbulence generated at a single source, 
the representative length scale may be relatively smal! due to the transfer of energy 
from large to small eddies via cascading. However, as the turbulence decays, the 
energy transfer also decreases. This causes the intensity of the smaller ed@icsuaae 
decrease faster than the larger eddies. Thus, despite the formation of smaller eddies by 
cascading. larger eddies predominate; as the smaller eddies disappear due to viscosity, 
the dissipation rate decreases and the representative length scale increases. [Ref. 10: p. 
48] 


1S 


2. Boundary Layer Effects 

In view of the fact that small-scale turbulent length scales can affect the 
missile boundary layer development. the effects of such interactions should be 
reviewed. An important factor to be considered is the point at which asymmetric 
vortices begin to form. More significant, however, is the point at which any type of 
flow separation occurs, be it sv~mmetric or asymmetric. Delaying all types of flow 
separation will inevitably decrease the tendency toward the formation of asymmetric 
vortices. 

A dimensionless parameter useful in analyzing the boundary layer effects 1s the 
skin friction coefficient, C,. In general. a fluid which 1s moving over a relatively 
smooth surface (low C,) will experience flow separation sooner than would the same 
flow moving over a rough surface. This is a result of the turbulent flow close to the 
rough surface. Such a boundary layer will resist separation longer due to the high 
kinetic energy which reduces the tendency toward flow deceleration at that point. 
Therefore, increasing C, generally delays the onset of flow separation. 

Ee MIMEMiamcata presented by tiancock and Bradshaw (see Reference 11) 
and Meier and Kreplin (see Reference 12) show an inverse relationship between 
Pamoulence length scales and the skin friction coefficient: as the length scale is 
decreased, C, increases. More specificaliy, the maximum value of the skin friction 
coefficient occurs when the length scale is on the order of the boundary laver thickness 
6. Both of these experiments used turbulence grids in wind tunnels, which ideally 
generate isotropic homogeneous turbulence at a large distance from the plane of the 
grid [Ref. 10: p. 48]. 

Thus, as the ratio of dissipation length scale to boundary layer thickness 
(L,’0) approaches unity, the onset of flow separation is delayed. In other words, the 
presence of small-scale turbulence (of a dimension comparable to the boundary layer) 


decreases the tendency to form asymmetric vortices. 


3. Marine Environment 
The atmospheric boundary layer into which the vertically launched surface-to- 
mmenussiic (VLSAM) is launched is best described as turbulent [Ref. 13: p. Sl}. In 
general, the turbulence fluctuations are greater over land than over sea. However, the 
characteristics of turbulence in the marine environment are not fully understood, and 


the effects on missile boundary laver development deserve investigation. 


iy 


The atmospheric boundary laver (ABL) is formed by the interaction of the 
atmosphere and the surface (land or sea) over which it flows. Turbulence in this region 
is caused by a transfer of heat, momentum, and miass. The lowest segment of the ABL., 
termed the surface layer, 1s characterized by mechanically produced turbulence 
resulting from surface roughness. The surface laver comprises the lower ten percent of 
the ABL, and is on the order of 50 meters in height. Furthermore, tne majority oe 
flow in the surface laver itself can be considered horizontally homogeneous. 
{Ref. 1: pp. 4-5] 

A measure of the general roughness of the surface 1s given by Z°uiie 
roughness length. Its vaiue is determined as a function of the mean wind Velocity at 
various elevations above the surface of interest. Combining the roughness length with 
the elevation and windspeed, both the turbulence intensity and length scale can be 
determined empirically [Ref. 14: p. 5]. Typical surface roughness lengths for the ocean 
are 0.001 < Zs 0.01, where Ze is measured in meters; at a ten meter elevation, this 
can result in a turbulence intensity range between 13 and 17 percent (see Reference 14 
ele 

Apparently, the turbulence intensity in the marine surface laver can be 
significant. However, the effect of such turbulence fluctuations is highly dependent 
upon the length scales present. For the range of roughness lengths (0.001 < z < 
0.01) and the elevation (z = 10 m) given above, the length scale 1s approximately 85 
meters (see Reference 14 p. 8). This value represents a scale very much larcersinamed 
conventional missile, and as such would have little effect on its boundary laver 
development. Yet. the cascade effect presented above allows for a decrease™imem. 
length scale until viscous forces dissipate the energy. Therefore, it is possible that 
length scales which are initially much larger (for example, approximately 85 meters) 
than the dimension of a missile would decrease, or “cascade,” to scales where thev 


could affect the development of the missile boundary laver. 
p ry 1ay 


4. Mach and Reynolds Number Effects 
The effect of the Mach number M and the Reynolds number Re on the 
development of the missile boundary layer and asymmetric vortices 1s dependent largely 
on the geometry of the body. Exhaustive research (see References 4. 7, 8, and 6) has 
been conducted, often with mixed results. The many different combinations of slender 
noses and afterbodies make comparisons of data very difficult. However, some 


generalities are possiole. 


With regard to the Mach number, the magnitude of the side force Y tends to 
decrease as M is increased. The degree to which tne force is reduced, as well as the 
point at which the reduction becomes significant, is determined by the shape of the 
body. At one extreme, a blunt-nosed cylinder experiences a rapid drop in ¥ as the 
flow approaches transonic velocities. This is caused by the tendency to form a 
separation bubbie at the nose of the blunt cylinder (see Figure 1.1). [Ref. 4: p. 99] 

On the other hand, a tangent ogive-nosed cylinder may undergo asymmetric 
vortex shedding into the supersonic region. Due to the lack of a separation bubole, 
strong side forces develop which can affect the body significantly. In both the blunt- 
and ogive-nosed cases, however, the side force Y tends to decrease with increasing 
iiememuniber [Ref. 7: p. 5]. 

The effect of the Reynolds number is best summarized by Ericsson and Reding 
[Ref. 4: p. 107] ; the vortex-induced side force reaches a maximum at a critical 
Reynolds number, when conditions are fully subcritical on one side of the body and 
fully supercritical on the other. A word of caution is necessary; the results for both the 
Mach and Reynolds number effects are based on several different experiments. 
Consistencies in data which hold for several body configurations would seem to 
indicate a valid trend. However, slight changes in test conditions often reverse the 
previously observed trends. Therefore, the tendency to draw broad conclusions in this 
area can be severely misleading. Przirembel and Shereda present an enlightening view 


of this problem in Reference 15. 


C. EFFECTS ON VERTICALLY-LAUNCHED SURFACE-TO-AIR MISSILES 

The use of vertically-launched surface-to-air missiles (VLSAM) is currently 
receiving wide acceptance. The advantages presented above demonstrate distinct 
improvements in reliability and flexibility. However, the flow regime into which a 
VLSAM is launched is not fully understood. As a result. the aerodynamic effects are 
difficult to predict. The effect of small-scale turbulence on asymmetric vortices has not 
yet been exaniuned experimentally. Research into this area is the goal of this thesis. 

It 1s usually desired that experimental research have general application in an 
area of current interest. For this reason, a generic VLSAM model, similar to current 
missile designs now in use, was developed by the author for this thesis. This generic 


design will hereafter be referred to as “the VLSAM” or “the model.” 


1. Launch and Crosswind Velocities 
In order that such research prove to be of current interest, is is necessary to 
examine the relative velocities of the VLSAM at launch and the crosswise 
experiences. Although each type of missile has a different launch velocity vs. time 
profile. a generic launch profile is used here to demonstrate relative velocities.” [ie 


generic data presented in Table i were provided by Farley (see Reference 1G)iaiieae 


data in Table 1 are presented graphically in Figure 1.4. 
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assumes 10 g missile longitudinal acceleration (generic) 


Examination of Figure 1.4 vields a slope of +322 ft'sec*. This corresponds to 
the acceleration in the vertical direction, a,. Substitution into kinematic equations at 
time t = 0.2 seconds, with a launcher height Zz, = 30 ft above the ocean surface, yields 
a height z ~ 56 ft. Thus, while still in the surface layer (Zz ~ 30 m), the VESA Vili 
a velocity of V, = 164 ft/sec. 

The assignment of a crosswind velocity is somewhat more arbitrary. It Is 
dependent on the ambient wind speed, but also on the speed of the launch platform 
itself. In simplifving the determination of the crosswind speed, a maximum and 
ininimum speed can be selected. With the minimum obviously at zero, the maximum 1s 
simply the sum of the mean wind speed and the ship’s speed. For a mean windspeed of 
20 m’sec [Ref. 1: p. 64] and a ship speed of 20 knots, the crosswine speed” V Wages 
iicece 1 — oseit sce). 

[In a two-dimensional cartesian coordinate system, let V, = velocity in the x- 
direction, and V, = velocity in the z-direction. Vector addition of the components 


vields a resultant speed V. = 191 ft/sec at an angle @ of 31° from the z-direction (see 
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iemre 1.5). Ihe direction represents the angle of attack a@ which the VLSAM 
experiences while still in the marine surface layer only 0.2 seconds after launch. This 
places the VLSAM in Regime III, the asymmetric vortex region, almost immediately 


miter clearims the exit plane of the launcher. 


2. Other Launch Considerations 
There are many factors which affect the aerodynamics of a missile during the 


launch pnase of its flight. Some are inherent to the design of the missile, while the 
effect of other factors is determined by the ship’s orientation at launch. In the former 
Gavecory are plume (or jet) effects of the missiles engine, blast effects of the vented 
exhaust gases, and activation of the flight control systems. The latter category includes 


shipboard roll, pitch, and yaw, and ship airwake turbulence. 
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Plume and blast effects can affect the VLSAM in two wavs. The frsmuismmme 
direct effect of the propulsion svstem gases or exhaust gases on the aerodynamic 
surfaces of the missile. Additionally, the exhaust gases can directly impact the flowfield 
into which the missile is launched, especially if the gases are vented upward intomeme 
vicinity of the accelerating mussile. 

The manner in which the control systems respond to changes in mussile 

rientation 1s another factor to be considered. Various svstems are emploved; some 
use a simple ballistic flight during the boost phase, while others attempt to maintain 
the missile’s orientation to a reference plane. Matters are further complicated by 
accounting for the motion of the ship as the missile travels along the guide rail prior to 
clearing the launcher. 

In addition to the turbulence generated in the marine atmospheric boundary 
layer, turbulence is also present in the airwake generated bv flow over the ship's 
superstructure. The effects of this turbulent ship airwake on the VLSAM are largely 


unknown. Full understanding of the effects of the airwake on the missile boundary 


layer cannot be accomplished, however, without first studying the complex nature of 
the airwake itself. Research by Healey (see Reference 14) outlines the need for 
comprehensive experimental work which combines the oscillatory motion of the 
launching platform with the bluff body turbulence generated by the superstructure. 

The above considerations, while not forming a complete description, 
emphasize the multiple factors which impact upon the flight behavior of a missile 
during the launch phase. Except for the limited discussion in the preceding paragraphs, 


these additional effects wiil not be included in this thesis. 


HW. EQUIPMENT AND PROCEDURES 


A. APPARATUS 

The major pieces of equipment utilized in the experimental portion of this thesis 
consisted of a wind tunnel, a missile model with associated support structures, 
turbulence-generating grids, and data acquisition hardware software. Information 
concerning the specifications, construction, and configuration of all equipment used is 
provided in this section. 

1. Wind Tunnel 

The wind tunnel utilized for the experimental data was a low-speed, single- 
return horizontal tunnel installed in Halligan Hall at the Naval Postgraduate School in 
Monterey, California (see Figure 2.1). It 1s powered by a 100 horsepower elegime 
motor coupled to a three-blade variable pitch fan via a four-speed Dodge truck 
transmission. Turning vanes are installed at all four corners, with two small-mesh wire 
screens installed upstream of the settling chamber. A heavy wire screen is also installed 
downstream of the test section to prevent damage to the fan blades or turning vanes in 
the event of model failure during operation. [Ref. 17: pp. 3-1 to 3-7] 

The test section of the tunnel measures 45 inches by 32 inches, with corner 
fillets which house the test section lighting; the cross sectional area 1s 9.88 square feet. 
The contraction ratio of settling chamber area to test section area 1s approximately ten. 
The walls of the test section diverge slightly to allow for boundary laver growth 
without a reduction in the freestream pressure along the test section. As the test 
section is designed to operate at atmospheric pressure, a circumferential breather slot is 
installed downstream of the test section to replenish air lost through leaks in the tunnel 
walls. The tunnel was designed to provide velocities up to 290 feet per second in the 
fest section. ([Ret, 17; pps 5-4) 00) O20) 

A reflection plane is installed in the floor of the test section, which decreases 
the vertical dimension to 28 inches. In the center of the reflection plane is a flush- 
mounted turntable for adjustments in pitch angle. The angle of the turntable 1s 


remotely controlled with an electric motor installed beneath the tunnel. 
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Figure 2.1 Schematic of Wind Tunnel. 


2. Turbulence Grids 

Four turbulence grids were used in the tunnel. Ideallv, such grids would 
generate isotropic homogeneous turbulence over a range of scales and intensities. All 
four grids were constructed by the author, and were designed for rapid instaliation and 
removal to facilitate maximum data acquisition with minimum tunnel occupation time. 
Each grid was mounted in a wood frame, which was placed 73 inches ahead of the 
pitch axis of the model support system (see Figures 2.2 and 2.3). The grids were 
installed in the tunnel using flush-mounted bolts, which extended through the frame 
and tunnel walls, and were securely fastened from outside the tunnel. The 
specifications of the four grids are listed in Table 2. 

The wooden turbulence grids used were square-mesh, square-bar biplane grids 
(see Figure 2.4), while the wire grid was also a square mesh but with round bars. 


Experience reveals the biplane grid to be superior for the generation of nearlv-isotropic 
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TAB GE 
GRID SPECIPIGATAGs 


Grid M * da Md Material 
AO Se eee 
l a0 ac 5 wood 
2 SS i 5 wood 
3 2 0.5 5 wood 

4 Kae 0.0625 16 Wire 


* M = mesh width 


** d = bar diameter 
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Figure 2.2 Position of Turbulence Grids in Wind Tunnel. 


homogeneous turbulence. Both biplane grids with rectangular bars, and monoplane 
grids, have been shown to produce nonuniform and unsteady flow, possibly caused by 
the increase in the separated region downstream of each grid intersection [Ref. 11: p. 


285]. A drawback of grid-generated turbulence is the fact that it is only 
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Figure 2.3 Turbulence-Gencrating Grid Mounted tin Wind Tunnel. 


Mears. - isotropic. Castro [Ref. 18: p. 300] shows the longitudinal rms velocity 
ccmponent to be up to five percent greater than the lateral components. Tunnel 


Prockace with the wooden grids was 20 percent, and 6.25 percent with the wire grid. 


3. Missile Model and Support System 

ence, es vemodelmwas cesigned bythe author to be representative of current 
vertically-iaunched surface-to-air missiles. [It is intended to model a cruciform tall- 
Memtro! mssiig with long dorsal fins (see Figure 2.5); however, all of the “control” 
Surfaces are fixed in position. It was constructed of 6061 and 2024 aluminum alloy bv 
Naval Postgraduate Scnool personnel. It 1s designed to permit force. moment. and 
pressure measurements while being operated in a subsonic wind tunnel. The force and 
moment measurements are accomplished using an internal strain-gage balance mounted 
On a sting. Pressure measurements are made using pressure taps iocated just aft of the 
Semee nose. Due to the dimensions of the model, pressure ard force moment 


Measurements cannot be made concurrently. The nose, tail, and dorsal fins are 
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Figure 2.4 Square-Viesh Turbulence-Generating Grid. 


removable, both for access to the measurement devices as well as for future 
configurational changes. The dimensions of the model are summarized below: 

*) “Diameter = 41.7 sciies 

e Length = 22.85 inches 

e Missile length’diameter ratio = 13.06 

®  Ogive nose length = 4.0 inches 

e Ogive length diameter ratio = 2.29 

eS WWoralsean = 25-4 ooineles 

® Tail span = 5.50 inches 

e Center of pressure location ~ 13.5 inches iioinese nie 

The support svstem for the VLSAM missile model was adapted to this wind 

unnel by the author; its design and construction allow testing of different models of 
various configurations. Rigidly mounted to the turntable and supported at the tunnel 


ceiling, it permits simultaneous variation of the roll and pitch (or vaw) angles. 
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+. Hot-Wire Anemometry Equipment 
Hot-wire anemometry was utilized to obtain data for longitudinally mapping 
the turbulence generated in the wind tunnel. Equipment was installed in the tunnel for 
both the not-wire calibration and the turbulence mapping. Extensive use was made of 


existing data acquistion equipment, including computer hardware and software. 


a. Cahbration Apparatus 
the hot-wire calibration svstem is diagrammed in Figure 2.6. Much of the 
equipment and computer software used in the calibration was developed for prior wind 
Meme, research. A Dantec single-sensor wire probe (tvpe 55PI1) was installed between 
Boe pressure probes, One measuring total pressure and the other measuring static 
Meeoure (See Figure 2.7). The hot-wire pressure probe was mounted at the 
downstream end of the wind tunnel test section. The wire sensor was connected to a 


reir A-100 hot-wire anemonieter, while the pressure probes were connected via 
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transducers to a Paroscientific digital pressure computer. The hot-wire anemometer 
and the pressure computer were connected to a Hewlett-Packard PC Instruments 
system. consisting of a relay multiplexer, digital multimeter, digital input/output device. 
and a system power unit. The instrument system interfaced with an IBM Personal 


Computer AT. 
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Figure 2.6 Block Diagram of Hot-Wire Calibration Apparatus. 


Another major component in the hot-wire calibration system was an 


Interactive computer program designed to facilitate data acquisition (see Reference 19). 
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eure 2.7 Villot- Wire Sensor and Pressure Probes. 


Minor modifications to this program were made bv the author. The program, CALIB1 
(as modified), 1s contained in Appendix A. 
b. Turbulence Mapping Apparatus 
The apparatus used for the actual turbulence mapping is diagrammed in 
Figure 2.8. {The hot-wire and anemometer were the same as used in the calibration. 
Additional equipment included a Tektronix 465M oscilloscope and a Fluke 8050A 
digital muitimeter for acquiring true rms values. The hot-wire probe (see Figures 2.9 
and 2.10) was mounted on an aluminum cvlinder, approximately centered horizontally 
mee eitically in the wind tunnel. The probe could be positioned between 5 inches and 
Momiiciies tipstream of the cylinder, which was rigidly mounted to the floor of the wind 
tunnel. Additionally. the cylinder itself could be positioned longitudinally in the wind 
tunnel. Thus, the hot-wire sensor could be positioned in a region extending from the 
downstream end of the test section to approximately six inches downstream of the 


turbulence grids. 
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Figure 2.8 Block Diagram of Turbulence Mapping Apparatus. 


B. PROCEDURE 
1. Hot-Wire Calibration 

The ourpose of the hot-wire calibration was to generate a curve for correlating 
voltage fluctuations to velocity fluctuations. The procedure used for hot-wire 
calibration closelv followed that outlined in the TSI anemometer instruction manual 
and the hot-wire experiment manual [Ref. 19]. The following is a summary of the 
procedure utilized. 

The data acquisition system: was configured as shown in Figure 2.6. The hot- 
wire; pressure probe was mounted on a ngid aluminum stand at the downstream end of 
the test section. A shorting probe was inserted in the probe support and cable 
resistance was measured and recorded bv the anemometer. 

The wire sensor probe was then installed in the probe support and probe 
resistance R,,, was measured. The operating resistance R of the probe was then 


calculated using Equation 2.1, 
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Figure 2.9 Hot-Wire Sensor Mounted For Turbulence Mapping. 


Pe Gang Rog le = ie) (eqn 2.1) 
where @,, and R,, are supplied with the wire sensor and I. 1s selected by the operator. 
The calculated value of R was entered into the hot-wire anemometer. The wind tunnel 
was started and placed in operation at a speed (100 knots) above that expected in the 
turbulence mapping. The anemometer Wheatstone bridge and cable compensation 
controls were adjusted, according to the instruction manual, for maximum frequency 
response. The hot-wire calibration program CALIBI was then initiated on the 
microcomputer. [he wind tunnel was set to various speeds, using the 
mucromanometer, which extended through the range to be used with turbulence grids in 
place. The following parameters were recorded for each speed set in the wind tunnel: 
anemometer output voltage, test section velocity, tunnel temperature, dvnamuic 
pressure, static pressure, and air density in the test section. In addition, hot-wire 


voltage and wind tunnel velocity were stored in a data file for use in generating the 
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Figure 2.10 Hot-Wire Sensor Downstream of Turbulence Grid. 


hot-wire calibration curve. This calibration procedure was conducted both prior to and 


after obtaining the turbulence mapping data. 


2. Turbulence Mapping 

The purpose of the turbulence mapping was to determine both the turbulence 
intensity and length scales as a function of distance downstream of the turbulence- 
generating grids. It was assumed that the biplanar, square-bar, square-mesh grids 
generated nearly-isotropic homogeneous turbulence at a large number of mesh widths 
from the plane of the grid. Therefore, single-wire hot-wire measurements were only 
made in the center of the tunnel cross-sectional area. 

“A hot-wire calibration was conducted prior to the commencement of the 
turbulence mapping. The hot-wire probe support and guide tube were mounted on the 
aluminum cvlinder and positioned in the tunnel. The data acquisition system was then 
configured as shown in Figure 2.8. The cable resistance previously determined was 


entered into the hot-wire anemometer. The wire sensor probe was inserted in the 


probe support, and the probe operating resistance R was calculated and entered into 
tie anemometer. 

The digital multimeter provided an rms value of voltage fluctuations across 
the wire sensor. As the value displayed changed five times every two seconds, an 
average based on the values displayed was calculated. Thus, the operator would record 
the values displayed, then compute the “average” rms voltage fluctuation depending on 
the number of values recorded. The output voltage was recorded from the display on 
the hot-wire anemometer, and the signal was monitored on the oscilloscope display. 

An initial run was made with no grid in place. The wire sensor was at 
Position 0 (see Table 3) and 100 rms voltage readings were taken. Subsequent rms 
voltage averages were made with 50 values. This represented an observation period of 


approximately 20 seconds. 
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The following turbulence mapping runs were conducted with the grids in 
position. Runs were made for each grid with the wire sensor at each of the positions 
listed in Table 3. Once again, averages were computed for the rms voltage 
fluctuations, while the output voltage was obtained directly from the anemometer 


displav. During each run, the micromanometer setting and the tunnel temperature 


were recorded. The atmospheric pressure was checked and recorded periodically. A 
hot-wire calibration was conducted following the completion of the turbulence 


mapping. 


fil. RESULTS 


A. HOT-WIRE ANEMOMETRY 

The results of the hot-wire anemometry are separated into two areas; hot-wire 
calibration and turbulence mapping. 

1. Hot-Wire Calibration 

Hot-wire calibration was conducted both before and after the turbulence 
mapping. The results of the two calibrations, presented in Table 4 and Figure 3.1, 
show a slight divergence between the initial and final runs. Linear equations were 
generated for each curve in the regions utilized in the turbulence mapping. For a given 
voltage, the initial and final equations differed by approximately seven percent. The 
turbulence mapping was evenly distnbuted over a three day period; therefore, no 
inference could be made as to which was the “better” curve. Since rms values are 
presented as a fraction of the local mean velocity, a small error in the calibration has 
little effect. For this reason, an average of the two equations was used for all further 
calculations. 

In addition, it was determined that the remote thermometer, which was used 
as an input to the calibration computer program, was unrelabie. As a result, a direct 
reading thermometer was used to correct the velocity values generated by the 
calibration program. The corrections are reflected in Table 4 and Figure 3.1. 

2. Turbulence Mapping 

The turbulence mapping yielded a large amount of data which was used in the 
determination of turbulence intensities and length scales. As discussed in Chapter I], 
the rms voltage resulting from velocity fluctuations was averaged from a large number 
of readings (at least 50 for each run), each reading itself being a mean of values 
performed internal to the true rms multimeter. A statistical analvsis of the rms 
readings at Position 0 yielded a mean standard deviation of less than 6 percent of the 
average rms value for each grid. The hot-wire anemometry system electrical noise was 
less than 0.5 millivolts; as a result, its effect was negligible. The results of the 


turbulence mapping are presented in Tables 5, 6, 7 and 8, as well as in Figure 3.2. 
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TAB ieee 
HOT-WIRE CALIBRATION RES: 


PU EE FINAL 
Wire Voltage Velocity Wire Voltage Velocity 
Oe eee (see) __ _ (oa ig 

0.968 0.00 1.394 Gone 
1.278 46.64 1.426 73.59 
[57 Gie7e 1.447 80.62 
1.384 C25 1.470 87.14 
1.414 isan ot 98.81 
1.433 80.79 S27 LOCH 
1.452 85.80 1.545 116.44 
1.483 100.44 ce [2oes3 
1.510 108.94 [oe 132.98 
i323 116.14 1.592 142.97 
1.542 Lee 1.604 149.61 
Led 134.69 Lets Rink 
ono 141.68 L620 163.05 
1.581 147.43 1.629 170.71 
1.593 [56.52 1.634 L 732602 
1.60] 162.63 

1.607 169.34 

ol 175.24 

1.618 181.70 

1.624 190.16 
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Figure 3.1 Hot-Wire Calibration Results 


The turdulence of the tunnel, with no grid in place, was calculated with a 
mean voltage of 1.582 volts and an rms voltage of 0.86 mullivolts: the turbulence 
intensity of the tunnel was 0.23 percent. Each of the runs with the wood grids (Grids 
#1, #2, and #3) was conducted with a dynamic pressure q = 10.00 cm H,O. This 
resulted in a fairly consistent mean voltage. (The dynamic pressure was used only to 
ensure simular wind tunne! velocities for each run; with the grids in place the actual 
value of the dvnamuic pressure, determined from the static rings, was in error due to the 
grid frame disturbance.) The wire grid (Grid #4) had to be run at a significantly lower 
dynamic pressure, q = 5.45 cm H,O, to obtain mean voltages which were similar to 
those obtained with the wood grids. The difference was caused by the much smaller 
pressure loss across the wire grid (blockage = 6.25%) as compared to the wood grids 
(blockage = 20°%%). 
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TAREE 
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Sensor ox Mean Voltage RMS Voltage 
Oe ee (VOLS) gull 

8 6.88 uel 6055 

i 13.88 1.495 45.51 

6 20.88 ae 33.01 

> Zoul5 1.498 27.86 

4 52a 1.497 23.03 

3 bee 1.495 19.67 

2 50.63 Lon 15.41 

| 7 00 ae [372 

0 64.63 1.516 V2ES2 

A als 1.536 [oo 

B Su Eeee 10.70 

TABLE 6 
TURBULENCE DAVTASFOKIGER!® 22 

Sensor ex Mean Voltage RMS Voltage 
CO ae a ae (01S is Deg EUVONGS 

8 6.88 1.438 64.84 

7 13.88 1.469 39.66 

6 20.88 1.490 2a 

5 Zoele 1.490 PANIES: 

4 Sale 1.498 ee 

3 3713 1.508 13.36 

2 50.63 5 12 

| 57.63 Pos 11.44 

0 64.63 les 2 10.45 

A Foals bce 2 

B 82.13 12550 9.07 


TAI 7 
TURBULENCE DATA*FOR GRID #3 


Sensor Xx Mean Voltage RMS Voltage 
eee ees) J al Nokes) 

8 6.88 1.455 44.81 

7 13.88 1.489 25-00 

6 20.88 197 17 o7 

5 Uae eS 2 14.55 

4 322 oS ee 

3 39.13 1.510 10.28 

2 S008 1516 8.44 

l 573s 525 1205 

0 64.63 1.528 f30 

A fogs 12537 6.62 

B Sele 1.542 6.31 

TABLE §& 
Heo LeEaCeE DATA FORFGRID #4 

Sensor mx Mean Voltage RMS Voltage 
oom ee ee a ESOS) 

8 6.88 1.480 Soo 

i 13.88 1.490 107 

6 20.88 1.49] 4.30 

5 Poel 1.504 B-o2 

4 S203 1.504 3.06 

3 Seale 1.505 2G 

2 50.63 323 Deal 

57.63 1.524 1.94 

0 64.63 esl eo 

A Hoel 1.529 1.65 

B S23 oly oe 
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a. Turbulence Intensity 

The equation of the average calibration curve was used to determine the 
mean velocity of tne wind tunnel at each position and corresponding mean voltage 
listed in Tables 5 - 8. Additionally. the slope of the curve was calcuiated to yield the 
ratio of the change in velocity to the change in voltage. This ratio was calculated for 
each posiuon and rms voltage in the above tables. 

The hot-wire apparatus measured the resultant of the longitudinal and 
lateral components of the turbulence. For grid-generated turbulence, the ratio of the 
longitudinal rms intensity to that of the lateral component is approximately 1.05 
{[Ref. 18: p. 300]. This difference was accounted for in calculating the longitudinal 
turbulence intensity u’ from the measured velocity fluctuation P’ and the mean velocity 
L. Thus, the actual longitudinal turbulence intensity was obtained from the measured, 


two-component Velocity fluctuation using Equation 3.1 
u-U = 0.724 (PU) (eqn Sam) 


The calculated values for turbulence intensity are shown in Table 9 and Figure 3.2 (see 
also Appendix B). 

Inspection of Figure 3.2 shows somewhat erratic behavior close to the 
grids. Most notable is the magnitude of the turbulence intensity immediately 
downstream of Grids #1 and #2. The smaller-mesh Grid #2 generated higher if than 
the larger-mesh Grid #1. This discrepancy was most probably caused by the lateral 
position of the sensor with respect to the grids. For Grid #1, the sensor was positioned 
downstream of an “open” region between the square bars; for Grid #2, the sensor was 
downstream of a grid intersection. The data from Table 5 show a relatively high mean 
voltage for Grid #1 (Position 8), while Table 6 shows a relatively high rms voltage for 
Giid 22 at the same position: 

The above inconsistencies occurred within two mesh widths (2M) for both 
grids. Previous experimental data show a minimum of 6M 1s required for the flow to 
achieve reasonable homogeneity [Ref. 18: p. 300]. Thus, for all the grids evaluated, the 
data should not be considered homogeneous for the first 30 inches; this corresponds to 
Positions § through 4. The relative turbulence intensity of the four grids at Position 0, 
which corresponds to a position on the nose of the missile model, is shown in Figure 
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Sensor, iurpulence Intensity (percent) 

Position __ Grid#1___Grid#2__ Grid #3 _ Grid 
§ 17.42 20.66 Laas 1.80 
a [oe 12.08 aS 1.68 
6 peor 8.08 5.26 ze 
Ss 8.15 6.46 4.18 Os 
4 6.76 5026 ono 0.89 
3 a7 6 4.49 2.96 One7 
Z 4.44 3.60 2.40 0.60 
I aes S425 ae 55 
0 3.66 3.00 2.05 Qrot 
A aa Zero 1.84 0.46 
B oo) 2.54 1.74 0.45 


A final observation involves the “character” of the generated turbulence. 
For all data runs, with one exception, the oscilloscope displayed a “uniformly random” 
pattern of high-frequency turbulence. In other words, the voltage fluctuations 
appeared to be continuous and relatively uniform across the oscilloscope. However, 
for Grid #1, with the sensor close to the grid (Position 8), the voltage fluctuations were 
interspersed with periods of almost no disturbance. The appearance was that of high- 
frequency turbulent fluctuations punctuated by bursts of laminar flow, indicating 
“smooth” airflow between turbulent vortices. The disparity in turbulence between Grid 
#1 and Grid #2 is further explained by the appearance of these laminar bursts; as 
laminar flow increases, the rms voltage fluctuations decrease, which lowers the 


turbulence intensity. 


b. Turbulence Length Scale 
The turbulence length scales were calculated using the turbulence intensities 
at the recorded positions. Meier [Ref. 20: p. 7] described a method for correlating the 
decay of the turbulence intensity with the distance downstream of a generating grid. 
Furthermore, Hancock [Ref. 11: p. 285] and Castro [Ref. 18: p. 300] link the dissipation 


length scale L, with the decay of turbulence intensity. 
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Figure 3.2 Grid-Generated Turbulence Intensity Curves. 


: : . eee * 
Linear approximation of the curves ( T,,” ) 0-8 versus downstream distance 


“x” yields lines of the form ol eauationes, 
CT = Ae ae) (eqn 3.2) 


where A and B are constants to be determined for each grid [Ref. ll: p. 283)e 50s 


above linear approximations are combined with Equation 3.3 
‘ ES Pee z 
(du ed ae = eee (eqn 3.3) 


and solved simultaneously (Ref. 11: p. 285]. The solution of these equations yields 


Equation 3.4, 
Le/M = 0.8 (A769 yx M — By?” (eqn 3.4) 
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Figure 3.3. Turbulence Intensity at Position 0 in Wind Tunnel. 


relating the length scale and the longitudinal distance. Thus, the constants A and B 
Were obtained for each grid. Substitution into Equation 3.4 gave the length scale L,, 
presented in Table 10 and Figure 3.4 (see also Appendix B). 

These results show a gradual increase in the length scale with increasing 
distance from the grid. As presented in Chapter I, this increase is due primarily to the 
cascade effect. As the turbulent eddies “cascade” to smaller length scales, they decay 
More rapidly due te viscosity. In turn, the average’ length scale, in effect, becomes 
larger with increasing distance from the grid. Thus, the increase in length scales shown 
in Table 10 also indicates the presence of smaller, more rapidly dissipated length scales. 
It is the smaller length scales, on the order of the boundarv layer thickness, which 
affect the development of the boundary layer. The relative turbulence length scale of 
the four grids at Position 0, which corresponds to a position on the nose of the missile 


model, is shown in Figure 3.5. 
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TABI ae 
GRID-GENERATED LENGTHS @ ees 


Sensor Length Scale (inches) 

Posivon =, nda! One ee Sa eee 
§ LOZ 0.84 055 0.14 
/ LG Or7 0.64 0.16 
6 ee 1.06 Orne Outs 
5 1.34 LlZ 0.76 Ong 
4 1.43 120 G32 Oy 1 
3 ion 27 O37 On? 
2 1.64 1258 0.95 0.24 
l i 1o 1.44 0.99 0.25 
0 be Fal 1:03 0.26 
A 1.86 leo 1.09 O27 
B 1.91 1.62 a 0.28 


An additional observation from Table 10 and Figure 3.4 is the comparison 
between the bar diameter of the grid and the length scale at Position 8. At that point, 
all of the square-bar grids (Grids #1, #2, and #3) yield a length scale which 1s only ten 
percent greater than the respective bar diameter. These three grids all have a M,d ratio 
of five. On the other hand, the wire grid (Grid #4) 1s a round-bar grid with a Md ratio 


of sixteen; at Position 8, its length scale is more than twice the bar diameter. 
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Figure 3.4 Grid-Generated Length Scale Curves. 
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Figure 3.5 Turbulence Length Scale at Position 0 in Wind Tunnel. 





IV. CONCLUSIONS AND RECOMMENDATIONS 


This experimental study has established a basis for further research on a 
vertically-launched surface-to-air missile in a turbulent flowfield at high angles of 
attack. The development of a generic VLSAM model, as well as the data collected and 
analyzed during the wind tunnel turbulence mapping. serve as a foundation for the 
continuation of this experimental work at the Naval Postgraduate School. 

The turbulence mapping yielded several important points relevant to the effect of 
turbulence on the development of asymmetric vortices and the missile boundary layer. 
Of primary importance was the determination of the length scales for each of the grids. 
Calculations showed that the length scale increased with increasing distance from the 
grid. This was caused by the dissipation of the smaller length scales which result from 
the cascade effect. As the smaller eddies dissipate due to viscosity, the proportion of 
larger eddies becomes greater, effectively increasing the “average” length scale. 
fle wever, it is mot the increase wnich is of interest, but rather the difference in the 
magnitude of the scales generated by the various grids. Increasing the bar diameter d 
resulted in increased length scale L,. As the parameter of interest is the relation 
between the length scale and boundary laver thickness, such a variation would appear 
necessary to allow testing of the model under different flow conditions. 

The turbulence intensity was found to increase with increasing bar diameter. The 
intent is to model turbulence intensities which might realistically be encountered during 
the launch of a VLSAM. As with the length scales, variation of the turbulence 
intensities, between the extremes, 1s considered necessary for validating experimental 
results. In both instances ee and ie) however, it must be remembered that a suitable 
distance must be maintained from the grid to allow the flow to develop reasonable 
homogeneity; data collected near the generating grid may be affected by vortices 
trailing off of the grid bars. 

Although additional research and experimental studies of the VLSAM have 
already begun, recommendations for additional research are included. As noted above, 
the VLSAM model has been designed for force. moment and pressure measurements. 
Therefore, a full set of experimental studies should be conducted, to include 


force, moment data and pressure measurements at all roll angles and angles of attack. 


ol 


These experiments should be further supported by surface and flowfield flow 
visualization and hot-wire measurements of the vortex flowfield. This would allow 
correlation of the turbulence intensities and length scales with the boundary layer 
development, asymmetric vortex formation and the resultant side forces. Alterations 
and modifications should be made to the ogive nose and dorsal tail fins to examine the 
effect of such changes on the development. Finally, after obtaining extensive 
experimental data. the results should be compared with predictions or models of 


turbulent flowfields over slender bodies of revolution at high angles of attack. 


csr 
ho 


APPENDIX A 
HOT-WIRE CALIBRATION COMPUTER PROGRAMI 


The computer program CALIBI was utilzed in the format shown below. The 
program) shell, which precedes the main program and configures the instruments, 1s not 
listed. 

Peavestments tome Onicinal proeram (see Kelerence 19) were intended to 
eliminate options not required for the hot-wire calibration, as well as to improve the 
accuracy of the data obtained. The call to the subroutines for reading and selecting 
angle of attack (Lines 4940-5450) was deleted. Additionally, the arrays for the hot-wire 
voltage and pressure differential were enlarged to accomodate a greater number of data 
points for each reading. 

Aiter obtaining the calibration curves. however, 1t was discovered that several 
additional corrections were necessary. These included scaling the velocities obtained to 
reflect the use of a direct-reading thermometer. Additionallv, Line 4860 adjusts the 
dynamic pressure, obtained from the pressure probes mounted adjacent to the hot-wire 
sensor (see Figure 22/7), Using the correction factor (0.93) for the tunnel static rings. 
This adjustment was determined to be incorrect, and the data obtained from the 


program were corrected accordingly. 


1000 ' PROGRAM "CALIB1" 

1010 OPTION BASE 1 

1020 WANTAOA = 0 

CS: Omme 

LOso: 

1050 * SET HP BOXES TO INITIAL SETTINGS FOR DATARACOUTS IEE 

160" 

1070 

1080 FILES = "“HOTWIRE.HPC"” 

1090 CALL INITIALIZE-SYSTEM(FILES) 

1100 CALL ENABLE.SYSTEM 

1110 CALL SET.FUNCTION( DMM. 01, DECVOLTS) 

1120 CALL SET.RANGE (DMM.01, AUTOM) 

i130 CALL SET.SPEED( DMM. 01/7 R2. 5) 

1140 CALL ENABLE.INT. TRIGGER (DMM. 01) 

1141 CALL SET.NUM. BITS (DIG. IN2O1L, R16} 

1142 LEVEL = 5 

1143 CALL SET.THRESHOLD(DIG.IN.01, LEVEL) 

fio. 

11607 

TiO" PRINT] PERST SCREEN 

TSO" 

3 0. 

1200 BEEP 

lzZilQ CEs 

1220 COLOR 4: LOCATE 5,14: PRINT “4k ee Ak RRR KKRKKKK KKK KKK KKK KKK KKK KR RK RRR RK KK KK 
kkk RRM 

1230 LOCATE 9,14: PRINT "keke kkKKK KEK KKK KKK kk KKK KKK KKK KKK KK KEEN 
1240 COLOR 1: LOCATE 9,37: PRINT “PART 1": LOCATE 10,33: COLOR 4: PRINTS 33 


waa eK KN 


L250 
1260 
1279 
2S 
L290 
1309 


13.10 
1220 
1330 
1340 
1330 
1400 
1410 
1420 
1430 
1440 


LOCATE 6,14: PRINT “*": LOCATE 7,14: PRINT “*": LOCATE 8,14: PREND 
LOCATE 6,65: PRINT "“*":; LOCATE 7,65: PRINT “*™: LOCATE 8,652 (PREli a 
COLOR 1: LOCATE 7,18: PRINT "WELCOME TO THE HOT WIRE FLOW MEASUREMENT LAB" 
LOCATE 12,27: PRINT “THESIS BY LT GREG DOREMUS"™ 


LOCATE 15,20: PRINT "THIS PROGRAM IS USED IN CALIBRATING THE " 

LOCATE 17,13: PRINT "HOT WIRE SYSTEM FOR WIND TUNNEL AIR FLOW MEASUREMENTS" 
COLOR) 4: LOCATE 21,24 

INPUT "ENTER 1 TO CONTINUE, 2 TO ABORT";ANS1 

IF ANS1 <> 2 GOTO 1340 ELSE 5950 

CLS: LOCATE 10,36: COLOR 20: PRINT "WAIT...": COLOR 1 

1 

1 

' PRINT SECOND SCREEN 

1 

1 

CLS: LOCATE 6,10: PRINT "TO START THIS EXPERIMENT YOU MUST FIRST SET UP THE 


HOT WIRE" 

1450 LOCATE 8,12: PRINT "SYSTEM. THE HOT WIRE LABORATORY MANUAL WILL EXPLAIN THI 
Ss" 

1460 LOCATE 10,6: PRINT “PROCEDURE IN DETAIL. ONCE THE SYSTEM SET UP IS) 7GRaaas 

SATISFACTION" 

1470 LOCATE 12,10: PRINT "YOU MUST MANUALLY START THE WIND TUNNEL AND SET THE SL 
OWEST" 

1480 LOCATE 14,16: PRINT "TUNNEL VELOCITY YOU WILL USE IN YOUR EXPERIMENT.” 

1490 COLOR 4: LOCATE 20,24: PRINT “ENTER 1 TO CONTINUE, (20 1O scent 

1500 LOCATE 21,28: INPUT “OR 3 TO GO BACK ONE PAGE 7ANS2 

1510 COLOR 1: ON ANS2 GOTO 1560, 59501230 

£52 070% 

1530 * PRINT THIRD SCREEN —- START WIND fUNNEL 
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1540 


S50. 

Pete clLs: LOCATE 10,22: PRINT “IF HOTWIRE SYSTEM IS SET UP PROPERLY" 

1570 COLOR 4: LOCATE 12,29: PRINT "START WIND TUNNEL NOW" 

1580 LOCATE 20,24: PRINT "ENTER 1 TO CONTINUE, 2 TO ABORT" 

1590 LOCATE 21,28: INPUT "OR 3 TO GO BACK ONE PAGE"; ANS3 

1600 COLOR 1: ON ANS3 GOTO 1610,5950,1440 

1540 ' 

PSo50 * 

1560 CLS: LOCATE 10,22: PRINT "IF HOTWIRE SYSTEM IS SET UP PROPERLY" 

1570 COLOR 4: LOCATE 12,29: PRINT "START WIND TUNNEL NOW" 

PagOMLOCATE 20,24: PRINT “ENTER 1 TO CONTINUE, 2 TO ABORT" 

1590 LOCATE 21,28: INPUT "OR 3 TO GO BACK ONE PAGE"; ANS3 

1600 COLOR 1: ON ANS3 GOTO 1610,5950,1440 

1610 * 

620° ' 

1630 ' PRINT FOURTH SCREEN 

1640 ! 

re 50 ' 

1660 CLS: LOCATE 8,13: PRINT “A SERIES OF WIND VELOCITY MEASUREMENTS IS TO BE TA 
KEN" 

1670 LOCATE 10,12: PRINT "AT WIND SPEEDS VARYING FROM THE SLOWEST TUNNEL VELOCIT 
ve 

ieco LOCATE 12711: PRINT "TO BE USEDVIN Y@URPEXPERIMENT UP TO THE HIGHEST VELOCI 
yee 

1690 LOCATE 14,15: PRINT "A MINIMUM OF 10 READINGS SHOULD BE TAKEN TO CREATE" 
my OO LOCATE 16,29: PRINT “THE CALIBRATION CURVE." 

1710 COLOR 4: LOCATE 20,24: PRINT "ENTER 1 TO CONTINUE, 2 TO ABORT" 

1720 LOCATE 21,28: INPUT "OR 3 TO GO BACK ONE PAGE";ANS4 

#730 COLOR 1: ON ANS4 GOTOMwE790,5950, 1560 

1740 ! 

750° ' 

1760 ' ENTER DATA POINTS FOR CALIBRATION CURVE 

770 

so ' 

'790 CLS: SE@CATE 12,12: PRINDMYENTER THE YOU WILL TAKE TO 
CREATE” 

1800 COLOR 4: LOCATE 12,22: PRINT "NUMBER OF DATA POINTS": COLOR 1 


1810 
1820 
1830 
1840 
350 
(16) 
1860 
1870 
1880 
1890 
1900 
7210 
f220 
a 30 
1940 
£950 
1960 
1370 
1980 
no 90 
2000 
2010 


LOCATE 14,19: INPUT "THE CALIBRATION CURVE FOR THIS EXPERIMENT."; POINTS 

CLS 

TEST = 1 

CALL OUTPUT (RELAY .MUX.01, TEST) 

DIM Z(POINTS,2), X(POINTS), Y(POINTS), C(100), D(100), KNOTS(POINTS,2), DEC 


OPEN DATA FILES 


TIMES = "00:00:00" 

COLOR 4 

LOCATE 14,18: PRINT "CREATING DATA FILES FOR CALIBRATION VOLTAGES" 
COLOR 1 

OPEN “HOTWIRE.DAT" FOR OUTPUT AS #1 

CLOSE #1 

TIME = VAL(RIGHTS (TIMES, 2) ) 

IF TIME < 5 GOTO 1970 ELSE 1990 

$ 


' START DATA ACQUISITION 
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2020 5 

2030" 

2040 READING = 0: CHANNEL = 0 
2050 FOR READING = 1 TO POINTS 


2001 CLS: LOCATE 10,15: PRINT “DEPRESS seule): ON THE PARASCIENTIFIC TRANSDU 
CER" 

2032 COLOR 4: LOCATE 10,30: PRINT "Pi: 3ecorere 

20210 * 

2030." 


2040 READING = 0: CHANNEL = 0 

2050 FOR READING = 1 TO POINTS 

2051 CLS: LOCATE 10,15: PRINT “DEPRESS BUTTON ON THE PARASCIENTIFIC TRANSDU 
CER" 

Z20o2 COLOR 4;>" LOCATE 107,30; 9PRINT | Fl" Colona. 

2053 LOCATE 12,15: PRINT "MAKE SURE NO OTHER FUNCTION BUTTONS ARE DEPRESSES. 
2054 COLOR 4: LOCATE 21,28: PRINT “PRESS ENTER TO CONTINUE” 

zZ205> AS = INKEYS: IF A$ = "" THEN 2055 

2060 CLS: -LOCATE SS 717. COLOR 1 

2070 IF READING = 1 GOTO 2030 ELSE (2100 

2080 PRINT "IS THE PROBE POSITIONED FOR YOUR FIRST READING]. 

2090 GOTO 2110 

2100 LOCATE 8,19: PRINT "IS THE PROBE POSITIONED FOR A NEW READING?" 

2110 COLOR 4: LOCATE 9,36: INPUT "1 > YES” ANSS 

2120 IF ANS5S <> 1 GOTO 2060 

2130 GOLOR 1 

2140 CLS: LOCATE 16,19: COLOR 1: PRINT "IS THE TUNNEL SET FOR A NEW FLOW VELOC 
Tee 

2150 COLOR 4:° LOCATE 17,36: INPUT “i= VES Ansa 

2160 IF ANS6 <> 1 GOTO 2140 

2270 CLS: COLOR 20: LOCATE 6,35: "PRENT “Siam. 

2180 COLOR 1: LOCATE 13,24: PRINT "DATA VALUES FOR HOT WIRE VOLTAGE" 

2190 LOCATE 15,25: PRINT "AND TEST SECTION VEEGCT Vy eakre 

2200 LOCATE 17,30: PRINT "NOW BEING COLLECTED" 


25.16 FOR CHANNEL = 1 TO 2 

2220 ON CHANNEL GOTO 2230,2342 

2236 CALL OUTPUT (RELAY .MUX.01, CHANNEL) 
2240 CALL SET. RANGE (DMM. 01,R20) 

2250 I = 0: SUM1 = 0: FOR I = 1 TO 100 
2260 IF I = 1 GOTO 2270 ELSE 2290 
2276 TIME = .5 

2280 CALL DELAY (TIME) 

2290 CALL MEASURE(DMM.01,C(I)) 
2300 SUM1 = SUM1 + ABS(C(I)) 

2516 NEXT I 

2320 AVG1 = SUM1/100 

2330 Z(READING,CHANNEL) = AVG1 

2340 GOTO 2780 

2542 I=0 

2344 SUM2 = 0 

2 2c0 FOR I = 1 TO 100 

2410 CALL MEASURE(DIG.IN.01,D(I)) 
2420 SUM2 = SUM2 + D(I) 

2430 NEXT I 

2440 DN = SUM2/100 

2442 GOSUB 5480 

2450 Z (READING, CHANNEL) = PRESS 

2460 GOSUB 4270 

2470 GOSUB 4700 

2480 Z(READING, CHANNEL) = VEL 

2490 KNOTS (READING, CHANNEL) = KTS 


2300 
Z2oL6 
Ziae 0 
2530 
2540 
2000 
2500 
2510 
2520 
2530 
2540 
2950 
2560 
2570 
Zoe0 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
27.0 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 


2900 
22 1G 
2920 
2230 
2940 
2250 
22.60 
2270 
2980 
2290 
3000 
S010 
3020 


BreeermGils -  b©CATE 2,32: PRINT “CALIBRATION DATA” 

LOCATE 4,17: PRINT "HOT WIRE VOLTAGE" 

LOCATE 4,46: PRINT "WIND TUNNEL VELOCITY": COLOR 4 

LOCATE 6,19: PRINT USING "#.####";Z(READING,1): COLOR 1 
LOGATE 6,26-SPRINT "VOLTS": COLOR 4 

LOCATE 6,41: PRINT USING "###.###"7Z(READING,2): COLOR 1 
BEEP: CBS: LOCATE 2,32: PRINT “CALIBRATION DATA" 

LOCATE 4,17: PRINT “HOT WIRE VOLTAGE" 

LOCATE 4,46: PRINT “WIND TUNNEL VELOCITY": COLOR 4 

LOCATE 6,19: PRINT USING "#.####";Z(READING,1): COLOR 1 
LOGATE 6,26: PRINT "VOLTS": COLOR 4 

LOCATE 6,41: PRINT USING "###.###"7Z2(READING,2): COLOR 1 
LOCATE 6,49: PRINT "FT/SEC": COLOR 4 

LOCATE 6,58: PRINT USING "“###.###" ; KNOTS (READING,2): COLOR 1 
LOCATE 6,66: PRINT “KNOTS S. SEOCATE (6,56: PRINT “=" 

EOCATE 8,15: PRINT “THiS, IS) DATA, POINT “7; READING;" OF ™; POINTS 
LOCATE 10,32: PRINT "WIND TUNNEL DATA" 

LOCATE 12,15: PRINT "TUNNEL TEMP =" 

COLOR 4: LOCATE 12,29: PRINT USING “##.##";TF: COLOR 1 
LOCATE 3127.35: “PRINT. “DEG Fr” 

LOCATE T4715: PRINT “fUNNEL STATIC PRESSURE =" 

COLOR 4: LOCATE 14,40: PRINT USING "####.##";PSTAT: COLOR 1 
LOCATE 14,48: PRINT "LB7Fr- 2" 

LOCATE 16,15: PRINT "TUNNEL AIR DENSITY =" 

COLOR 4: LOCATE 16,36: PRINT USING "#.#####";DENSITY: COLOR 1 
LOCATE 16744: PRINT. "LB/ FTL" 3" 

LOGATE 28715: PRINT “DYNAMIC PRESSURE =" 

COLOR 4: LOCATE 18,34: PRINT USING "##.##";Q: COLOR 1 

LOCATE 18,40: PRINT “LB/FT“2" 

IF READING = POINTS GOTO 2750 

COLOR 20: LOCATE 21,24: PRINT "SET NEW FLOW VELOCITY IN TUNNEL" 
GOLOR 4: “LOCATE 23,27: PRINT "PRESS ANY KEY TO CONTINUE" 
Nom—sONKEYS: IF AS = "" THEN 2760 

COLOR 1 


NEXT CHANNEL 
NEXT READING 


BNUM1 
BNUM2 
BNUM3 
BNUM4 
1 


' 
' STORE 
' 
' 


Ze) 
Z (POINTS ,2) 
KNOTS (1,2) 


RNOTS (POINTS, 2) 


DATA IN DATA FILE LABELED "HOTWIRE. DAT" 


CLS: LOCATE 8,15: PRINT "DATA VALUES ARE BEING STORED IN DATA FILE LABELED" 


COLOR 4 

I = 0: 

FOR I = 
A(T) 
Y (1) 
Ci 


Oo 
} 
uo ue tl 


a 
tm 
rs 
ry 
Hi 


¢ 


e 
Cc 
D 
D 


LOCATE 10,34: PRINT “HOTWIRE.DAT": COLOR 1 
Lt = 0: C2aa— 0: DI = 0: D2°= 0 
eho ONS 

SOR(Z(I,2)) 

(750s Ae) ae 
1 + X(T) 
2ervx( be 2 
Lote Cla 
Zee Lee (1) 


OPEN "HOTWIRE.DAT" FOR APPEND AS #1 
FOR J = 1 TO POINTS 
PRINT #1,USING"###4.###8"7Y¥(J), X(J) 


wo =  O: 


ae 


SOSO0ONEXT oJ 

3040 CLOSE #1 

3050) ' 

3060 ' 

3070 ' SOLVE FOR SLOPE AND Y-INTERCEPT OF REGRESSION LINE 
3080 '! 

3090 ! 

3030 NEXT @ 

S0405CLOSH et. 

3050" 

3060 
3070 
3080 
3090 
3100 
2110 
2120 
S130 
3140 
3150 
3160 
3170 
3180 
3190 LOCATE 16,12: PRINT "WOULD YOU LIKE TO SEE A SUMMARY OF YOUR CALIBRATION DA 
TAM 

3200 LOCATE 18,33: PRINT "ON THE SCREEN?": COLOR 4 

3210 LOCATE 20,27: INPUT "ENTER 1 FOR YES, O FOR NO";ANS7 

3220 IF ANS7 <> 0 GOTO 3230 ELSE 3640 

3230 COLOR 1: CLS: LOCATE 2,32: PRINT "CALIBRATION DATA" 

3240 LOCATE 4,9: PRINT "HOT WIRE VOLTAGE (VOLTS)" 

3250 LOCATE 4,41: PRINT "WIND TUNNEL VELOCITY (FT/SEC)": COLOR 4 

3260 J =o 

3270 FOR J = 1 TO POINTS 

3280 LOCATE J+5,18: PRINT USING"#. ####":Z(J,1) 

3290 LOCATE J+5,52: PRINT USING" ###.###"2Z(7,2) 

3300 NEXT J 

a3 10° 

Sa 210! 

3330 ' SEND CALIBRATION DATA TO HP LASER JET PRINTER 

8346.5! 

3350 * 

3360 COLOR 1: LOCATE 22,20: PRINT "WOULD YOU LIKE A HARD COPY OF THIS DATA?" 
3370 COLOR 4: LOCATE 23,27: INPUT "ENTER 1 FOR YES, 2 FOR NO";ANS8 

3380 IF ANS8 <> 2 GOTO 3390 ELSE 3640 

3390 CLS: LOCATE 12,26: PRINT "DATA SENT TO LASER PRINTER": COLOR 1 

3400 TIMES = "00:00:00" 

3410 TIME = VAL(RIGHTS (TIMES, 2) ) 

3420 IF TIME < 5 GOTO 3410 ELSE 3430 

3430 LPRINT: LPRINT 

3440 LPRINT TAB(32): LPRINT "CALIBRATION DATA" 

3450 LPRINT 


SOLVE FOR SLOPE AND Y-INTERCEPT OF REGRESSION LINE 


= (POINTS*D2-C1*D1) / (POINTS*C2-C1*2) 
F ABS(A) < 1E-08 THEN A=0 

= (D1-A*C1)/POINTS 

F ABS(B) < 1E-08 THEN B=0 


PRINT SUMMARY OF CALIBRATION DATA ON SCREEN 


-*-= = = «~HMWHPY = - = «= 


3460 LPRINT TAB(9): LPRINT "HOTWIRE VOLTAGE (VOLTS) WIND TUNNEL VELO 
CITY" 

3470 LPRINT TAB(42): LPRINT "(FT/SEC) (KTS)" 

3480 LPRINT 

3490 J = 0 


3500 FOR J = 1 TO POINTS 

S510 LPRINT TAB(18); 

3520 LPRINT USING"#.####";7Z2(0,1); 
3530 LPRINT TAB(43); 
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3540 LPRINT USING" ###.###";2(3,2); 

3550  LPRINT TAB(61); 

3560 LPRINT USING'###.###";KNOTS(J,2) 

3570 NEXT J 

3580 LPRINT CHRS$(12) 

3590 ' 

3540  LPRINT USING" ###.###";2(J,2) 3 

3550 LPRINT TAB(61); 

3560  LPRINT USING"###.###"; KNOTS (J,2) 

3570 NEXT J 

3580 LPRINT CHR$(12) 

3590 ! 

3600 ! 

3610 ' PRINT VIEW PLOT OF CALIBRATION CURVE ON THE SCREEN 

3620 ! 

3630 ! 

3640 COLOR 1: CLS: LOCATE 12,13: PRINT "WOULD YOU LIKE TO SEE A PLOT OF THE CALI 
BRATION CURVE?" 

3650 COLOR 4: LOCATE 14,27: INPUT "ENTER 1 FOR YES, 2 FOR NO";ANS9 

3660 COLOR 1: IF ANS9 <> 2 GOTO 3670 ELSE 3730 

3670 CLS: SHELL "COMMAND/C HOT1.BAT" 

3680 ' 

3690 ' 

3700 ' SEND HARD COPY OF CALIBRATION CURVE TO HP LASER JET PRINTER 

a0 ' 

3720 ' 

3730 COLOR 1: CLS: LOCATE 12,19: PRINT "WOULD YOU LIKE A HARD COPY OF THIS GRAPH 
ot 

3740 LOCATE 14,27: COLOR 4: INPUT "ENTER 1 FOR YES, 2 FOR NO";ANS10 

3750 COLOR 1: IF ANS10 <> 2 GOTO 3760 ELSE 3840 

3760 CLS: COLOR 1: LOCATE 10,21: PRINT "THIS PROGRAM WILL AUTOMATICALLY LEAVE" 
3770 LOCATE 12,22: PRINT "BASICA AND ENTER A GRAPHICS PROGRAM." 

3780 LOCATE 14,14: PRINT "THE PLOT TAKES APPROXIMATELY 4-6 MINUTES TO COMPLETE" 
3790 COLOR 4: LOCATE 21,27: PRINT "PRESS ANY KEY TO CONTINUE" 

3800 A$ = INKEY$: IF A$ = "" THEN 3800 

3810 COLOR 1 

3820 SHELL "HOT2.BAT" 

3830 BEEP 

3840 OPEN "CHAIN.DAT" FOR OUTPUT AS #1 

3850 PRINT #1,USING"####.###"; A, B, BNUM1, BNUM2, BNUM3, BNUM4, TF 

3860 CLOSE #1 

3870 ' 

3880 ' 

3890 ' PRINT CONCLUDING SCREEN AND SUMMARY OF CALIBRATION RANGE 

3900 ' 

3910 ' 

3920 CLS: LOCATE 8,14: PRINT "THIS CONCLUDES THE CALIBRATION CURVE PORTION OF YO 
UR" 

3930 LOCATE 10,10: PRINT "EXPERIMENT. THE HOT WIRE SYSTEM IS NOW ACCURATELY CALI 
BRATED" 

3940 LOCATE 12,13: PRINT "WITHIN A WIND TUNNEL VELOCITY RANGE OF" 

3950 COLOR 4: LOCATE 12,52: PRINT USING" ###.###";Z(1,2) 

3960 LOCATE 12,61: PRINT "FT/SEC": COLOR 1 

3970 LOCATE 14,17: PRINT "(": LOCATE 14,18: PRINT USING "###.###";KNOTS (1,2) 
3980 LOCATE 14,26: PRINT "KTS) TO": COLOR 4 

3990 LOCATE 14,34: PRINT USING "###.###";Z(POINTS,2): LOCATE 14,42: PRINT "FT/SE 
Cc": COLOR 1: LOCATE 14,49: PRINT "(" 

4000 LOCATE 14,50: PRINT USING "###.###";KNOTS (POINTS, 2) 

4010 LOCATE 14,58: PRINT "KTS)" 

4020 COLOR 4: LOCATE 20,24: PRINT "ENTER 1 TO CONTINUE, 2 TO ABORT" 


=), 


4030 LOCATE 21,28: INPUT "OR 3 TO GO BACK ONE PAGE“ VANS 

4040 ON ANS11 GOTO 4100,5950,3730 

4050 ! 

4060 ' 

4070 ' PRINT LAST SCREEN AND CALIBRATION EQUATION 

4080 '! 

4030 LOCATE 21,28: INPUT "OR 3 TO GO BACK ONE PAGE";ANS11 

4040 ON ANS11 GOTO 4100,5950,3730 

4050 ! 

4060 ' 

4070 ' PRINT LAST SCREEN AND CALIBRATION EQUATION 

4080 ! 

2090! 

4100 CLS: LOCATE s6)30: PRINT “Yo=(" 

4110 LOCATE 8,34: PRINT USING "“##.###";A: LOCATE 8,40: PRINT "xX" 

4120 IF B>=0 THEN LOCATE 8,42: PRINT “+": LOCATE 8,44: PRINT USING “##.###";B EL 
SE LOCATE 8,42: PRINT "=": LOCATE 8,44: PRINT USING "###.###" 7;ABS(B) 

4130 COLOR 1: LOCATE 10,10: PRINT "THIS IS THE STRAIGHT LINE EQUATION OF STR Eae2 


IBRATION CURVE” 


4140 
4150 
4160 
4170 
4180 
459 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 


LOCATE 12,19: PRINT “IT WILL AUTOMATICALLY BE ENTERED IN PART 2" 
LOCATE 14,26: PRINT “OF YOUR HOT WIRE EXPERIMENT" 

COLOR 4: LOCATE 20,27: PRINT “ ENTER 1 TO START PART TWO" 

LOCATE 21532: PRINT "ENTER 2 TO TABOR 

LOCATE 22,25: INPUT “OR ENTER 3 TO GO BACK ONE PAGE. Atlee 

COEOR 

ON ANS12 GOTO 4210,5950,3920 

LOAD "PART2",R 

' 


' SUBROUTINE FOR MEASURING WIND TUNNEL TEMPERATURE 


Al = 38.709457# 
A2 = .037085566# 
A3 = 5.649552E-05 
BO = .10086091# 
Bl = 25727.94369# 
B2 = -767345.8295# 
B3 = 78025595.81# 
B4 = -9247486589# 
B5 = 6.97688E+11 
B6 = ~2.66192E+13 
B7 = 3.94078E+14 
B8 = 0 

B9 = 0 


REF.OUT = 8 
CALL OUTPUT (RELAY.MUX.01,REF.OUT) 
CALL ENABLE.OUTPUT (RELAY.MUX.01) 
CALL SET. FUNCTION (DMM.01,DCVOLTS) 
CALL SET.RANGE(DMM.01,R2) 
CALL DISABLE.INT.TRIGGER(DMM.01) 
COUNT =\.5 
TIMES = "00:00:00" 
IF TIMER < COUNT GOTO 4480 ELSE 4490 
CALL MEASURE(DMM.01,V) 
TR = 100*V 
ER = (TR*(A1+TR* (A2+TR*A3)))*10*~6 
CALL SET.RANGE(DMM.01,R200MILLI) 
TC =a 
CALL OUTPUT (RELAY.MUX.01,TC) 


60 


4550 
4560 
4570 
4580 
4590 
4600 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4722 


COUNT = .5 

TIMES = "00:00:00" 

IF TIMER < COUNT GOTO 4570 ELSE 4580 
CALL MEASURE (DMM.01,ET) 
E = ER + ET 
Z = BS+E*(B6+E* (B7+E* (B8+E*B9) ) ) 

COUNT = .5 

TIMES = "00:00:00" 

IF TIMER < COUNT GOTO 4570 ELSE 4580 
CALL MEASURE (DMM.01,ET) 
E = ER + ET 
Z = BS+E* (B6+E* (B7+E* (B8+E*B9) ) ) 
T = BO+E*(B1+E* (B2+E* (B3+E* (B4+E*Z) ))) 
TF = (1.8*T) + 32 
RETURN 


' 
t 
' SUBROUTINE FOR WIND TUNNEL DATA: 

: STA@EC PRESS IN TEs? secrion 

: AIK DENSI TyYeIN TEST SECTION 

; DYNAMIC PRESSURE IN TEST SECTION 

: TEST SECTION VELOCITY, FT/SEC AND KTS 
t 


CLS: COLOR 1: LOCATE 10,15: PRINT "DEPRESS BUTTON ON THE PARASCIENTIFIC 


TRANS DUCER" 


4724 
4726 
4728 
4730 
4800 
4802 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 


COLOR 4: LOCATE 10,30: PRINT "P2": COLOR 1 

LOCATE 12,15: PRINT "MAKE SURE NO OTHER FUNCTION BUTTONS ARE DEPRESSED" 
COLOR 4: LOCATE 21,28: PRINT "PRESS ENTER TO CONTINUE" 
AS = INKEYS: IF AS = "" THEN 4730 

CALL MEASURE (DIG.IN.01, DN) 

GOSUB 5480 

PSTAT = PRESS*144 

TTEMP = TF + 459.67 

GASCONST = 53.3 

DENSITY = PSTAT/(GASCONST*TTEMP) 

DELTAP = Z(READING, CHANNEL) *144 

Q = DELTAP/.93 

IF Q < 1.5 THEN Q = 0 

Gc = 3.174 

VEL = SQR(ABS((2*Q*GC) /DENSITY) ) 

KTS = VEL * .5924 

RETURN 

t 


t 
' SUBROUTINE FOR READING AOA 
t 
t 


INPUT.CH = 5 

CALL OUTPUT (RELAY .MUX.01, INPUT.CH) 
CALL ENABLE.OUTPUT (RELAY .MUX.01) 
TIME = 3 

CALL DELAY (TIME) 

CALL MEASURE (DMM.01,A0A) 

HAVEAOA = (AOA-.003) *1000/10 
RETURN 

t 


1 
' SUBROUTINE FOR AOA SELECTION 
t 
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5090 


5100 WHILE (CINT(ABS(WANTAOA - HAVEAOA) * 10) / 10) >= 
ABSDIFF = ABS(WANTAOA ~ HAVEAOA) 
IF ABSDIFF > 36 THEN ABSDIFF = 36 
SWITCH.1 = 1 


S10 
5220 
5130 
5090 
5100 
S10 
52.0 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
3210 
5220 
5230 
5240 
5250 
5260 
S270 
5280 
5290 
5300 
5s 10 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
Bozo 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 


WHILE (CINT(ABS(WANTAOA - HAVEAOA) * 10) / 10) >= 
ABSDIFF = ABS(WANTAOA - HAVEAOQA) 
IF ABSDIFF > 36 THEN ABSDIFF = 36 
SWITCH.1 = 1 
SWITCH.2 = 2 
CALL OPEN. CHANNEL (RELAY.ACT.01,SWITCH.1) 
CALL OPEN. CHANNEL (RELAY.ACT.01,SWITCH. 2) 
IF WANTAOA > HAVEAOA GOTO 5210 ELSE 5340 


IF (CINT(ABSDIFF*10)/10) >= .4 THEN TIMING = 


ao 


225 


THIS PART OF THE SUBROUTINE MOVES THE TURNTABLE FORWARY 


COUNT = ((ABSDIFF * TIMING) + 190) / 1000 
IF COUNT > 3 THEN COUNT = 3 

CALL CLOSE. CHANNEL (RELAY.ACT.01,SWITCH. 1) 
CALL DELAY (COUNT) 

CALL OPEN. CHANNEL (RELAY.ACT.01, SWITCH. 1) 


TIME = 45 
CALL DELAY (TIME) 


GOSUB 4940 
GOTO 5440 


WEND 


RETURN 


OUNT 


DEC 


7-“-*- 7 ©-FQ = + -= = 


DEC (16 
DEG (ie 


DEC(6) 
DEC (2) 
t 

IF DN 
COUNT 


100 ELSE TIMING 


" THIS PART OF THE SUBROUTINE MOVES THE TURNTABLE IN REVERSE 


IF (CINT(ABSDIFF*10)/10) >= .4 THEN TIMING 


COUNT = ((ABSDIFF* TIMING) + 190)/1000 
IF COUNT > 3 THEN COUNT = 3 

CALL CLOSE.CHANNEL(RELAY.ACT.01,SWITCH.2) 
CALL DELAY (COUNT) 

CALL OPEN. CHANNEL(RELAY.ACT.01,SWITCH. 2) 


LIME. = 2.5 

CALL DELAY (TIME) 
GOSUB 4940 

GOTO 5440 


= 0: PRESS = 0 


DN IS THE MEASURED DECIMAL NUMBER 


IS THE DECIMAL TRANSFORM MATRIX 


PRESS IS THE CONVERTED NUMERIC PRESSURE VALUE 


) = 10: DEC(15) = 8: DEC(14) = 4: DEC(13) 

) = .8: DEC(10) = .4: DEC(9) = .2: DEC(8) 
= .04: DEC(5) = .02: DEC(4) = .01: DEC(3) 
= .004: DEC(1) = .002 


< 0 THEN PRESS = DEC(16) ELSE GOTO s5640 
= -32768! 
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100 ELSE TIMING 


De GDEC (12 = 
-1: DEC@in— 
8.000001E-03 


a 


08 


1 


— 
= 


SUBROUTINE FOR CONVERTING DECIMAL PRESSURE VALUES TO NUMERIC VALUES 


al 


5640 
5650 
5660 
5670 
5680 
5690 
5700 
5710 
97:20 
5730 
5740 
57.50 
5760 
5770 
5780 
5720 
5800 
5810 
5820 
5830 
5840 
SaO0 
5860 
5870 
5880 
5890 
5900 
52 LO 
ao 2.0 
Boo. 
5940 
S950 
5960 


IF (DN-COUNT)/16384 >= 1 THEN PRESS = PRESS + DEC(15) ELSE GOTO 5660 
COUNT = COUNT + 16384 
IF (DN-COUNT)/8192 >= 1 THEN PRESS = PRESS + DEC(14) ELSE GOTO 5680 


COUNT = COUNT + 8192 


IF (DN-COUNT)/4096 >= 1 THEN PRESS = PRESS + DEC(13) ELSE GOTO 5700 


COUNT = COUNT + 4096 


IF (DN-COUNT)/2048 >= 1 THEN PRESS = PRESS + DEC(12) ELSE GOTO 5720 


COUNT = COUNT + 2048 


IF (DN-COUNT)/1024 >= 1 THEN PRESS = PRESS + DEC(11) ELSE GOTO 5740 


COUNT = COUNT + 1024 
IF (DN-COUNT) /512 >= 
COUNT = COUNT + 512 
IF (DN-COUNT)/256 >= 
COUNT = COUNT + 256 
IF (DN-COUNT) /128 >= 
COUNT = COUNT + 128 
IF (DN-COUNT)/64 >= 
COUNT = COUNT + 64 
IF (DN-COUNT) /32 >= 
COUNT = COUNT + 32 
IF (DN-COUNT)/16 >= 
COUNT = COUNT + 16 
IF (DN-COUNT)/8 >= 1 
COUNT = COUNT + 8 

IF (DN-COUNT)/4 >= 1 


COUNT = COUNT + 4 
IF (DN-COUNT)/2 >= 1 
COUNT = COUNT + 2 

IF (DN-COUNT) = 1 
RETURN 

STOP 

CLS: SYSTEM 

END 


1 THEN PRESS = PRESS + DEC(10) ELSE GOTO 5760 
1 THEN PRESS = PRESS + DEC(9) ELSE GOTO 5780 
1 THEN PRESS = PRESS + DEC(8) ELSE GOTO 5800 
1 THEN PRESS = PRESS + DEC(7) ELSE GOTO 5820 
1 THEN PRESS = PRESS + DEC(6) ELSE GOTO 5840 


PRESS + DEC(5) ELSE GOTO 5860 


1 THEN PRESS 


THEN PRESS = PRESS + DEC(4) ELSE GOTO 5880 
THEN PRESS = PRESS + DEC(3) ELSE GOTO 5900 
THEN PRESS = PRESS + DEC(2) ELSE GOTO 5920 


THEN PRESS = PRESS + DEC(1) 
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APPENDIX B 
TURBULENCE MAPPING RESULTS 


The combined results of the turbulence mapping were presented in Chapter III. 
Shown below are the individual graphs which represent the turbulence intensity and 
length scale data for the four turbulence-generating grids. Attention should be given to 


the scales of the graphs. 
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TURBULENCE INTENSITY (PERCENT) 


LENGTH SCALE (INCHES) 


4 6 § 10 12 14 16 18 


eS 


0.4 0.8 ee 1.6 


0.0 





_#A 


GRID #1 


20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 


= 
Ss 


GRID #1 


0.0 20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 


Figure B.1 Turbulence Intensities and Length Scales For Grid #1. 
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TURBULENCE INTENSITY (PERCENT) 


LENGTH SCALE (INCHES) 


ne 


10 12 14 16 18 20 


0.8 1.2 1.6 ZO 2 4 6 8 


0.4 


0.0 





GRID #2 





0.0 20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 


re eee 


GRID #2 


0.0 20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 


Figure B.2 Turbulence Intensities and Length Scales For Grid #2. 
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TURBULENCE INTENSITY (PERCENT) 


LENGTH SCALE (INCHES) 


10 AZ 14 


8 


6 


& 4 


0.4 0.8 lee 1.6 2.0 0 


0.0 


GRID #3 
=) 
0.0 20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 
GRID #3 
0.0 20.0 40.0 60.0 80.0 


DISTANCE FROM GRID (INCHES) 


Figure B.3 Turbulence Intensities and Length Scales For Grid #3. 
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0.0 
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1.0 


GRID #4 


0.8 


0.6 


0.4 


0.2 


ee 


0.0 20.0 40.0 60.0 80.0 
DISTANCE FROM GRID (INCHES) 


LENGTH SCALE (INCHES) 
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Figure B.4 Turbulence Intensities and Length Scales For Grid #4. 
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